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Aba  tract 

The  objective  of  this  thesis  Is  to  minimize  the  terminal  position 
error  of  an  entry  vehicle.  To  do  this,  the  concept  of  minimizing  the 
position  sensitivity  coefficients  is  employed.  With  this  in  mind,  two 
purposes  are  established:  first,  to  investigate  a  set  of  optimal  solu¬ 
tions  minimizing  a  criterion  function  of  sensitivity  coefficients;  and, 
second,  to  form  a  suboptimal  guidance  law  for  these  trajectories, 
the  trajectories  are  generated  by  using  a  type  of  control  which  in¬ 
cludes  parameters  describing  the  vehicle  configuration  and  attitude. 
This  is  called  a  lumped  parameter  control.  Due  to  the  use  of  a  lumped 
parameter  control,  a  Valentine's  procedure  is  employed  to  restrict  the 
range  to  a  set  of  realistic  design  values.  A  closed-loop  optimal 
implicit  guidance  scheme  including  a  State  enHtnatrrr  -fa  n«i»H  -tn  a  a-lnm- 
latlon  of  selected  trajectories  to  test  the  validity  of  the  feedback 
gains  determined  for  the  trajectories. 

The  minimum  sensitivity  coefficient  trajectory  appears  to  be  one 
in  which  the  vehicle  first  levels  out  to  nearly  horizontal  flight  and 
then  dives  to  achieve  the  desired  terminal  flight  path  angle.  The 
values  of  the  radius  and  flight  path  angle  sensitivity  coefficients  for 
this  type  of  trajectory  are  about  an  order  of  magnitude  lower  than  for 
other  types  considered,  while  the  maximum  velocity  sensitivity  coef¬ 
ficient  remains  about  the  same. 

For  most  of  the  cases  considered,  there  appears  to  be  a  definite 
relationship  between  the  sensitivity  coefficients  and  range  angle  and 
rate  of  change  of  flight  path  angle.  For  all  cases  but  one,  the  flight 
path  angle  sensitivity  coefficient  is  nearly  a  linear  function  of 
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rug*  angle. 

Th*  simulation  r**ults  using  a  sat  of  feedback  gains  ud  assuming 
noisy  observations  drives  the  system  to  within  S  a,  2  ape,  and  0.2  deg 
at  the  terminal  point  with  Initial  errors  of  5000  m,  100  mps,  and  1.0 
deg  in  position,  velocity,  and  flight  path  ugle  respectively. 
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SYNTHESIS  OF  A  SUB OPTIMAL  GUIDANCE  LAM  TO  MINIMIZE 
THE  TERMINAL  ERROR  OF  AN  ENTRY  VEHICLE 

I.  TntrnrinrHftti 

Background 

When  a  vehicle  enters  the  earth's  atmosphere,  any  error  in  its 
position  vector  or  velocity  vector  will  cause  a  subsequent  position 
error  at  the  desired  terminal  point.  The  quantities  which  describe 
the  sensitivity  of  the  position  error  at  the  terminal  point  are  called 
sensitivity  coefficients  or  influence  coefficients.  One  method  of 
minimizing  the  position  error  at  the  terminal  point  is  to  minimize 
these  sensltlblty  coefficients.  A  trajectory  which  minimised  the  sens¬ 
itivity  coefficients  is  also  the  trajectory  which  provides  the  minimum 
position  error  at  the  terminal  point.  This  concept  was  first  consider¬ 
ed  by  Johnson  in  his  dissertation  (Ref.  1). 

Purpose 

The  basic  objectives  of  this  thesis  are  as  follows: 

A.  To  synthesize  and  analyze  a  set  of  optimal  open-loop 
trajectories  for  an  entry  problem  which  will  minimize  a  function  of  the 
sensitivity  coefficients  described  above. 

B.  To  generate  a  set  of  suboptimal  feedback  gains  which 
will  drive  the  vehicle  to  the  desired  terminal  point  in  the  presence 
of  small  errors  in  the  vehicle  position  and  velocity  vectors  at  any 
time  during  entry. 

In  theory,  the  use  of  a  set  of  feedback  gains  would  appear  to 
invalidate  thi  objective  of  minimizing  the  sensitivity  coefficients, 
since  with  perfect  measurements  and  an  optimal  feedback  scheme,  the 
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terminal  conditions  would  bs  analytically  satisfied.  However ,  in  prac¬ 
tice,  neither  of  the  two  assumptions  are  valid.  First,  consider  the 

environment ;  the  measurements  ere  extracted  rrtno  an  Inertial  mwun- 
ment  Unit,  which  is  in  turn  adjusted  by  the  overall  attitude  of  the 
vehicle.  Errors  from  the  environment,  then,  can  be  additive,  and,  at 
present,  must  be  compensated  for.  Conalder,  alao,  the  feedback  scheme; 
the  ideal  controller  would  require  aome  nonlinear  relationship  between 
the  control  vector  and  the  position  and  velocity  vectors.  Since,  in 
an  entry,  control  must  be  applied  in  near  real  time,  the  nonlinear 
relationship  is  impractical.  Thus  a  linear  relationship,  valid  for  a 
small  region  about  a  nominal  trajectory,  la  used,  and  the  errors  In¬ 
herent  in  an  approximation  are,  again  significant.  In  practice,  then, 
there  is  ample  motivation  to  seek  a  trajectory  where  the  influence  of 


I  i  .1 
atiuiibUt 


The  open- loop  controller  which  achieves  the  desired  optimal  tra¬ 
jectory  is  defined  to  include  all  possible  parameters  of  the  non-thrus¬ 
ting  entry  vehicle  controller  and  all  the  parameters  which  describe 
the  vehicle  configuration.  This  is  done  for  the  purpose  of  allowing 
the  vehicle  designer  maximum  flexibility  in  choosing  a  combination  of 
parameters  which  meet  the  control  requirements  and  which  also  fall 
within  reasonable  design  limits.  This  type  of  control  is  called 
lumped  parameter  control. 


Assumptions 

It  is  assumed  that  the  entry  control  vector  synthesized  yields  a 
planar  trajectory  with  respect  to  a  spherical,  non-rotating  earth. 
Since  the  control  used  here  Is  a  lumped  parameter,  control,  it  is  nec- 
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assary  to  restrict  th«  rang*  of  control  to  raaliatic  design  limit*. 

Th«  specific  limit*  used  herein  are  given  in  Chapter  IZZ. 

in  conacructiu*  u’r.i  feedback  echeme.  it  is  assumed  that 

«  linear  relation  between  the  control  vector  and  the  vectors  for  posi¬ 
tion  and  velocity  can  be  used.  It  is  further  assumed  that  the  use  of 
the  position  and  velocity  vectors  alone  in  the  feedback  scheme,  with  no 
use  being  made  of  the  sensitivity  coefficients,  is  sufficient  to  drive 
the  vehicle  to  the  desired  terminal  conditions.  It  is  necessary  to 
make  this  assumption  because  the  errors  in  the  sensitivity  coefficients 
cannot  be  computed  onboard  the  vehicle. 

Approach  To  The  Problem 

The  optimal  open-loop  trajectory  and  the  control  required  are 
found  by  first  variation  optimization  techniques  baaed  on  the  calculus 
of  variations.  Using  this  trajectory  and  control  as  a  set  of  nominal 
values  a  set  of  linear  perturbation  equations  are  solved  to  refine  the 
values  such  that  the  terminal  conditions  are  perclsely  satisfied. 

These  linear  perturbation  equations  are  solved  using  a  second  variation 
optimization  technique  based  on  the  calculus  of  variations. 

Oilce  an  open-loop  trajectory  and  control  are  generated  that 
satisfy  the  terminal  conditions,  a  set  of  linear  perturbation  equations 
are  solved  to  produce  a  set  of  feedback  gains  for  the  trajectory.  In 
order  to  observe  the  effects  of  implementing  the  feedback  control  in  a 
practical  situation  with  measurement  errors  and  noise  present,  a  com¬ 
puter  simulation  is  accomplished  using  the  previously  determined 
nominal  optimal  control  and  the  feedback  gains.  A  Kalman  filter  is 
therefore  applied  to  the  nonlinear  trajectory  state  equations  of  motion 
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to  estimate  tho  position  snd  velocity  vsctors  in  ths  presence  of  noisa. 

Ihs  technique's  used  in  both  the  open-  end  closed-loop  optimisation 
pewblvM  ptwmubvu  iu  utiv*  iii  Cll«|l  W<l£  Hi  Ala  U  XU  Ullp  (Vr  XX  4U1 
outline  for  the  Implementation  of  the  Kalman  filter,  and  a  method  for 
analytically  generating  the  initial  error  covariance  matrix  for  the 
filter  are  presented  and  discussed.  The  model  equations  to  be  used 
for  the  problem  of  minimizing  sensitivity  coefficients,  end  a  trans¬ 
formation  of  independent  variables  is  accomplished  in  Chapter  IX?.,  In 
Chapter  III,  the  derivation  of  the  lumped  parameter  control  is  also 
presented  in  depth.  The  application  of  the  techniques  given  in  Chapter 
II  to  the  equations  from  Chapter  III  is  found  in  Chapters  IV  and  VI. 

In  Chapter  IV,  the  open-loop  equations  are  derived  using  the  optimal 
theory  from  Chapter  II,  and  some  computation  difficulties  occurring 

during  flip  unpUrnMnn  nF  rtnHwl  nri»n-  1r»nn  nl  gnrH  Him*  ar»  i)4 «maaail . 

In  Chapter  VI,  a  closed-loop  simulation  algorithm,  implementing  the 
theory  for  the  Kalman  filter  is  presented.  Also  in  Chapter  VI,  a  deri¬ 
vation  for  the  closed-loop  feedback  scheme  is  presented.  The  results 
of  the  open-loop  part  of  this  thesis  are  presented  and  observations 
made  In  Chapter  V,  while  the  results  for  the  closed-loop  part  are 
presented  and  discussed  in  Chapter  VIII,  and  recommendations  concerning 
this  problem  as  well  as  areas  of  further  study  are  given  in  Chapter  IX. 
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II.  Optimisation  Te chnlguts 

The  purpose  of  this  chapter  is  to  present  a  summary  of  general 
optimal  control  techniques  and  of  optimal  state  estimation  techniques. 
The  runtrysgin  maximum  principle  and  the  necessary  conditions  for  opti¬ 
mality  associated  with  it  will  be  discussed  as  well  as  numerical 
methods  for  satisfying  these  necessary  conditions  to  determine  an  open- 
loop  control.  The  general  equations  for  an  optimal  linear  feedback 
control  will  be  derived  and  the  limitations  of  this  control  will  be 
discussed.  The  equations  for  a  Kalman  state  estimator  applied  to  a 
nonlinear  system  will  be  summarised  and  a  general  computational  algo¬ 
rithm  presented. 

Statement  Of  The  General  Optimal  Control  Problem 

The  mathematical  model  describing  the  dynamics  of  a  controllable 
process  can  generally  be  expressed  as  a  set  of  first-order,  nonlinear, 
vector  differential  equations  as  follows: 

X  -  f(x,u,t)  (2-1) 

where  the  vector  x  includes  all  the  quantities  which  are  necessary  to 
describe  the  dynamics  of  the  system  (e.g. ,  position,  velocity,  and 
flight  path  angle  for  an  entry  vehicle) ,  and  the  vector  u  Includes  all 
the  system  control  variables  (e.g.,  lift  control  and  drag  control  for 
an  entry  vehicle).  The  vector  x  is  referred  to  as  a  state  vector.  The 
initial  conditions  on  the  state  vector  are  denoted  as  follows: 

2£(t0)Ai°  (2—2) 
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The  final  conditions  ate  denoted  byt 

x(tf)  A.x£  (2-3) 

The  problem  under  consideration  is  that  of  choosing  a  control, 
from  some  acceptable  class  of  controls,  which  best  satisfies  some  sel¬ 
ected  criterion.  The  criterion  selected  is  in  general  to  minimise  some 
function  of  the  state  variables  and  the  control.  This  function  will  be 
referred  to  as  the  criterion  function  and  will  be  defined  as  a  scalar. 
The  control  which  minimizes  the  criterion  is  the  optimal  control  and 
is  denoted  u*.  The  selection  of  the  criterion  function  is  a  vital 
part  of  the  problem  since  it  specifies  the  desired  performance  of  the 
system.  The  general  criterion  function  treated  here  is  called  the 
Lagrange  form  and  is  expressed  as  fallows: 

J  -  /^f|t(x(T),u(T),T]dT  (2^4) 

The  next  section  discusses  the  conditions  under  which  J  is  a  minimum 
(or  a  maximum).  The  variable  $  is  some  function  of  the  states  and 
controls.  The  derivation  of  the  conditions  is  descuased  in  additional 
detail  in  references  10  and  11. 

Summary  Of  The  Pontryagin  Maximum  Principle 

The  function  $ls  assumed  to  be  continuous  and  differentiable 
through  the  second  order  with  respect  to  the  components  of  the  state 
vector  .x  and  with  respect  to  the  components  of  the  control  vector  u. 

By  using  Lagrange  multipliers  as  adjoint  variables,  the  criterion  func¬ 
tion  is  augmented  as  follows: 

J  “  /^f  {#[x(t)  ,u(x),x]  +Vr(T)^(x(x),u(T),T]  -  x(x)j}  dx  .(2-5) 
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It  should  be  noted  that  this  equation  la  basically  no  different  from 
the  original  criterion  function,  Eq(2-4).  Since  the  factor  multiply¬ 
ing  the  Lagrange  multiplier  in  *ig(£-5)  is  id"M  rally  equal  to  aero 
from  Eq(2-1). 

From  Eq(2-5),  a  scalar  function  called  the  Hamiltonian  ia  defined 
aa  follows: 

Hlx(t),u(t),Jk(t),t]  -  *[x(t),u(t),t]  +XT(t)f.[x(t)  ,u(t),tj  (2-6) 
Substituting  thia  into  Eq(2-5)  gives: 

J-  /tf(H(x(T),u(t),Xi(t),'t]  -  X,t(t)x(t))  dr  (2-7) 

fc0 

When  Eq(2-7)  la  integrated  by  parts,  the  result  ia: 

J  -  -  AT(T>X(T)|tf  +/tf{H(x(T),u(T),X(T),T]  -  AT(T)x(T))dT  (2-8) 

*0  S 

When  Eof2-8)  ia  minimised.  the  original  criterion  function  is 

minimised  and  the  state  vector  equality  constraints,  Eq(2-1)  are  Bat- 

* 

is fled.  In  order  to  minimise  Eq(2-8),  a  small  variation  is  made  in 
both  the  state  vector  and  the  control  vector  and  the  limit  is  taken  as 
the  small  variation  approaches  cero> 

When  the  small  variations  in  state  and  control  are  made,  the 
following  equation  is  obtained  (dropping  functional  notation): 

«J  -  t6xT(-A)]|tf  +/^f{6xT[|^  +  ll  +  6uT  ^  }dr  '  (2-9) 

—  t  co  ~  135  “  “  3U 

0 

When  the  limit  is  taken  as  the  small  variations  approach  zero, 
each  of  the  terms  in  Eq(2-9)  must  approach  zero.  This  gives  the  fol¬ 
lowing  necessary  conditions  for  optimality: 
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1 1  tf  "  0  (2~1C 

i-  -H 

X  -  -3-  -  f(x,u,t)  <2-13 

g  "  0  <2-12 

where  Eq(2-10)  is  the  transversallty  condition  (i.e. ,  boundary  condi- 

* 

tion) ,  Eq(2-ll)  is  the  set  of  adjoint  equations,  Eq(2-12)  is  the  ori¬ 
ginal  set  of  state  equations,  and  Eq(2-13)  is  the  gradient  equation. 
This  set  of  equations  constitutes  the  necessary  conditions  for  opti¬ 
mality,  but  they  are  not  sufficient.  Solution  of  these  equations 
leads  to  a  local  extremum  which  may  be  either  a  maximum  or  a  minimum 
but  cannot  be  considered  a  global  minimum. 

The  transversallty  condition,  Eq(2-10),  must  be  considered  for 
two  separate  cases.  When  the  initial  or  final  condition  on  any  com¬ 
ponent  of  the  state  vector  x  is  specified,  the  first  variation  of  the 
component  at  that  point  is  identically  zero .  When  the  Initial  or 
final  condition  on  any  component  of  the  state  vector  is  unspecified, 
the  corresponding  adjoint  variable  at  that  point  must  be  identically 
zero  and  the  first  variation  in  the  state  component  is  then  unconstr¬ 
ained.  This  is  illustrated  in  Table  I. 


Table  I.  Constraints  On  Variables 


State,  x 

Specified 

Unspecified 

Adjoint,  X 

Unconstrained 

0 

Variation,  <5x 

0 

Unconstrained 

GGC/EE/70-10 


Open-Loop  Optimal  Methods 

First  Variation  Gradient  Technique.  The  gradient  procedure  die- 
cussed  here  inherently  satisfies  the  state  and  adjoint  differential 
equations,  and  satisfies  the  gradient  equation,  Eq(2-13),  by  an  iter¬ 
ative  procedure.  The  procedure  also  inherently  satisfies  the  trane- 
versality  conditions.  From  Eq(2-9) ,  and  assuming  that  Eq(2-10) 
through  Eq(2-12)  are  satisfied,  the  following  relationship  may  be 
obtained: 

6 J  -  /"V  dT  (2-14) 

t0  ~ 

The  objective  of  this  procedure,  then,  Is  to  drive  the  gradient 
to  zero.  To  achieve  this  goal,  a  gradient  3H/du  is  computed  and  a 
change  In  control  6,u  calculated.  The  new  control  is  then  used  to 
compute  a  new  gradient  and  the  process  is  repeated  until  6u  approaches 
zero.  When  iu.  is  sufficiently  small,  the. gradient  is  near  zero  and 
the  last  necessary  condition  of  optimality  is  considered  satisfied. 

The  sign  of  the  required  change  in  control  can  be  found  from  the 
sign  of  the  gradient,  but  the  magnitude  is  more  difficult  to  determine. 
The  magnitude  is  generally  determined  by  searching  over  a  range  of  mag¬ 
nitudes  until  one  is  found  which  minimized  the  criterion  function. 

The  gradient  technique  is  derived  and  described  in  more  detail  in  ..ref 
erences  7  and  11 . 

Assuming  that  all  state  components  are  specified  at  the  initial 
time  and  unspecifigd^'at  the  final  time,  a  typical  computation  algori- 
thm  for  implementing  the  gradient  technique  is  as  follows: 

Guess  an  Initial  control  versus  time,  u. 

Integrate  the  state  equations  forward  from  the 


Step  1. 
Step  2. 

4T 
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given  initial  conditions  and  store  the  final  values. 

Step  3.  Using  zero  as  the  final  condition  on  the  adjoin ts 
and  the  final  conditions  on  the  states  from  step  2,  Integrate  the 
state  and  adjoint  equations  backward  and  store  the  computed  gradient 
versus  time.  The  sign  on  the  gradient  provides  the  sign  for  the  con¬ 
trol  change. 

Step  4.  Perform  a  one-dimensional  search  over  different 
magnitudes  of  6u  until  a  6u  is  found  which  minimizes  the  first  vari¬ 
ation  of  the  criterion  function  (dJ) . 

Step  5.  Construct  a  new  control  versus  time,  unaw  ■ 

Hold  + 

Step  6.  Co  back  to  step  2  with  the  new  control  values  and 
continue  until  the  gradient  is  sufficiently  small  or  until  the  change 
in  either  the  control  or  the  criterion  function  Is  sufficiently  smsll. 

First  Variation  Terminal  Error  Function  Technique.  The  gradient 
algorithm  described  above  always  satisfies  the  state,  adjoint,  and 
transversality  equations  and  works  to  satisfy  the  gradient  equation. 
The  terminal  error  function  technique  always  satisfies  the  state,  ad¬ 
joint,  and  gradient  equations  and  works  to  satisfy  the  transversality 
equation. 

At  each  end  of  the  problem,  there  are  2n  boundary  conditions 
where  n  is  the  number  of  states.  Of  these,  n  conditions  are  specifi¬ 
ed  and  n  are  unconstrained.  For  the  a  that  are  specified,  two  vectors 
are  formed,  z.  and  zjj,  where  z,  is  the  actual  value  of  the  states  or 
ad joints  at  the  boundary  and  ££  is  the  desired  or  specified  value.  A 
new  criterion  function,  Q,  is  defined  as  follows: 
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q  -  U  -  *o]Tl»  -  *d1  <2“15> 

where  Q  is  called  the  terminal  error  function.  When  this  function  is 
minimized,  the  transversality  conditions  are  considered  satisfied. 

Siuue  the  statu  ana  adjoint  equations  ana  the  gradient  equation  are 
always  satisfied,  every  solution  is  an  optimal  control  for  the  part” 
lcular  end  conditions  achieved. 

To  implement  this  technique,  the  states  and  adjoints  are  integrat¬ 
ed  simultaneously  from  either  end  of  the  problem.  The  integration  is 
performed  using  the  given  boundary  conditions  for  those  variables 
which  are  specified  and  estimated  or  guessed  values  for  those  varia¬ 
bles  at  the  other  end  of  the  problem,  the  estimated  values  are  changed 
and  the  process  is  repeated  until  the  integration  terminates  with  the 
desired  boundary  conditions.  ThiB  is  equivalent  to  transforming  the 
two-point-boundary-value  problem  (TPBVP)  into  either  an  initial  value 
problem  or  a  final  value  problem. 

A  necessary  condition  for  the  use  of  this  method  is  that  the 
control  be  completely  removed  from  the  problem.  For  this  to  be  done 
the  second  partial  of  the  Hamiltonian  with  respect  to  the  control  vec¬ 
tor  must  be  non-singular.  That  1b: 


If  the  determinant  is  not  zero,  the  control  can  be  expressed  as: 

u  "  h(XiA»fc)  (2-17) 

A  substitution  for  the  control  in  the  state  equations  can  be  accompli¬ 
shed. 

Once  the  unconstrained  boundary  conditions  are  estimated  and  the 
boundary  conditions  at  the  other  end  of  the  problem  are  found  to  be  in 
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arror,  a  gradient  la  computed  for  tha  criterion  function  with  respect 
to  the  estimated  boundary  conditions,  E^,  as  follows: 


VQ  .  Q(Ej  +  AEj)  -Q(Et) 
1  AE« 


(2-18] 


The  computation  of  thia  gradient  provides  a  search  direction  for  the 
change  in  boundary  conditions  which  will  satisfy  the  boundary  condi¬ 
tions  at  the  other  end. 

A  computational  algorithm  for  the  terminal  error  function  pro¬ 
cedure  is  as  follows: 

Step  1.  Construct  a  terminal  arror  criterion  function,  Q, 
corresponding  to  o&e  end  of  tha  problem. 

Step  2.  At  the  opposite  end  of  the  problem,  guees  the  n 
unspecified  values. 

3.  luLegiuie  both  che  scare  and  aajoinc  cuzrerentxax 
equations  using  the  specified  and  guessed  boundary  conditions  an  a 
starting  point. 

Step  4.  Determine  the  value  of  Q. 

Step  5.  Numerically  differentiate  Q  with  respect  to  the 

t 

guessed  boundary  conditions  to  form  V<2. 

Step  6.  Perform  a  one-demens ional  search  over  the  guessed 
boundary  conditions  in  the  direction  of  -7£  to  minimize  Q. 

Step  7.  With  new  guessed  boundary  conditions  from  Step  6, 
go  back  to  Step  3  and  continue  until  Q  is  sufficiently  small. 

The  only  difference  between  this  method  and  that  of  the  second 
variation  is  that  in  the  terminal  error  function,  a  linear  search  is 
performed  based  on  the  first  variation  equations;  whereas,  in  the 
second  variation  procedure,  numerical  perturbation  differential 
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equations  aze  solved  during  the  Integration  of  the  state  and  adjoint 
equations  to  determine  the  necessary  changes  In  the  guesses.  The  ter¬ 
minal  error  function  technique  described  here  is  due  to  Trushell  and 
Birta  given  in  Eefeneuue  8. 

Second  Variation  Technique.  In  the  open-loop  optimal  control 
technique  based  on  the  second  variation  of  the  Hamiltonian,  it  is  as¬ 
sumed  that  all  necessary  conditions  for  optimality  are  Identically 
,  satisfied  with  the  exception  of  the  boundary  conditions.  This  tech¬ 
nique  transforms  a  two-point-boundary-value  problem  (1. e. ,  a  problem 
In  which  some  conditions  are  known  and  some  are  unknown  at  each  end  of 
the  trajectory)  into  an  initial-value  or  final-value  problem  (i. e. ,  a 
problem  in  which  all  of  the  conditions  are  known  at  either  one  end  of 
the  trajectory  or  the  other). 

Since  the  gradient  is  defined  to  be  identically  aero,  it  is  as¬ 
sumed  here  that  the  control  vector  can  be  expressed  as  an  explicit 
analytical  function  of  the  state  and  adjoint  variables.  Based  on  this 
assumption,  the  control  function  can  be  substituted  into  the  state  and 
adjoint  differential  equations  to  completely  remove  control  from  the 
problem.  Then,  once  the  boundary  conditions  have  been  satisfied,  the 
final  open-loop  control  versus  time  can  be  computed  from  the  functional 
relationship. 

Based  on  these  assumptions,  the  state  and  adjoint  equations  are 
Integrated  starting  at  either  end  of  the  trajectory.  Errors  in  known 

values  at  the  end  indicate  the  initial  unknown  boundary  conditions  are 
in  error  and  must  be  changed.  The  entire  purpose  of  this  procedure  is 

to  determine  the  magnitude  of  these  unknown  boundary  conditions  such 
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that  tha  known  boundary  conditions  at  the  other  end  of  the  trajectory 
are  ■afclif'ted.  Thla  la  done  by  solving  a  set  of  perturbation  differ¬ 
ential  equations  simultaneously  with  the  state  and  adjoint  differential 
equations.  The  solution  to  this  set  of  perturbation  equations  provides 
the  sensitivity  of  errors  In  the  specified  values  at  one  end  of  the 
trajectory  to  changes  in  the  estimated  values  at  the  other  end.  These 
sensitivities  are  then  used  In  a  Newton- Raphson  procedure  to  change 
the  estimated  values. 

The  sensitivity  quantities  provided  by  this  method  will  be  refer¬ 
red  to  In  this  thesis  as  boundary  value  partlals  to  avoid  confusion 
with  the  position  sensitivity  coefficients  which  are  being  minimized. 
The  boundary  value  partlals  are  sensitivity  coefficients ,  but  the 
sensitivity  coefficients  being  minimized  relate  deviations  In  position 
■*1  ong  the  trajectory  tc  crrsrc  In  states  at  tha  tarsias!  paint,  wham 
as,  the  boundary  value  ratios  relate  errors  in  all  of  the  known  state 
and  adjoint  variables  at  one  end  of  the  trajectory  to  changes  In  the 
unconstrained  state  and  adjoint  variables  at  the  other  end. 

To  Illustrate  the  Newton-Raphson  procedure,  an  example  Is  consi¬ 
dered  in  which  the  first  q  states  of  a  total  of  n  states  are  rna-uned 

to  be  specified  at  the  Initial  point  and  all.n  states  etf  ^pd>‘.is'ij44 

! 

at  the  finAl  point.  Two  separate  cases  must  be  considered  to t  ;nla 
example.  One  case  deals  with  the  procedure  when  backward  integration 
is  selected  to  satisfy  the  initial  boundary  conditions,  the  other  with 
the  procedure  when  forward  integration  Is  selected  to  satisfy  the 
final  boundary  conditions. 

In  the  case  of  backward  integration  for  this  case,  the  unconstr¬ 
ained  adjolnts  at  the  final  point  are  Iteratively  adjusted  to  satisfy 
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the  specified  boundary  conditions  at  the  initial  point.  The  specified 
boundary  conditions  at  the  final  point  are  satisfied  because  these 
conditions  are  used  to  begin  the  backward  integration.  The  Newton- 
Raphson  equation  for  this  case  may  be  expressed  ns: 


6ju(to) 


*(0 


Out  C  t  o ) 

3X(tf) 


4X(tf) 


which  may  be  rewritten  as: 


SA(tf)  - 


3X(tf) 


-1 


34<t0) 

3X(tj) 

4^(t0) 

3y(t0) 

M  t0) 

(2-19) 


(2-20) 


where  ■  a  q-vector  of  the  devlaHnna  of  the  ntetee 

[jt( tQ )]  from  the  specified  values. 

6w(tQ)  ■  an  (n-q)  vector  of  the  deviations  of  the  initial  ad- 

joints  [A(tg)]  from  zero.  These  adjoints  are  specified 
to  be  zero  at  the  initial  point  because  the  correspond¬ 
ing  states  are  completely  unconstrained. 

6X^(tj)  -  the  necessary  change':  in  the  unconstrained  adjoint  var¬ 
iables  at  the  final  point  to  satisfy  the  initial 
boundary  conditions. 

Eq(2-20)  may  be  written  in  terms  of  an  iterative  equation  as  follows: 


A(tf)Q+1  -  A(tf)n  + 


3f(to^ 

-1 

0- 

3X(tf) 

(2-21) 
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where  <p°  is  Che  vector  of  specified  initial  states,  i(i(to)  is  the  vec¬ 
tor  of  actual  initial  states  from  the  backward  integration  and  oi(tg) 
is  the  vector  or  actual  Initial  adjolnts.  The  quantities  in  the  matrix 
to  be  inverted  may  be  derived  by  performing  a  Taylor  series  expansion 
on  three  of  the  necessary  equations  of  optimality.  These  equations 
are 


and 


i(£,u) 

(2-22) 

■  -  ih 

(2-23) 

dx 

-  0 

(2-24) 

The  partial  derivatives  of  these  three  equations  with  respect  to  A.(c j) 


are: 


ai(t) 

d 

M 

3x(t) 

■ 

3f(x,u) 

3A(tf) 

dt 

3A(tf) 

3A(tf) 

,02 

rt 

d 

3A(t) 

3 

3H 

[ 3A(tf) 

dt 

.3A,(tf) 

1 

,02 

rt 

3x(t) 

(2-23) 


(2-26) 


3 

3H 

A 

3G(t) 

dnd  3A(tf) 

3u(t) 

A 

3A(tf) 

(2-27) 


Performing  a  Taylor  series  expansion  on  these  three  equati  ns  gives: 


d_ 

3x(t) 

3f (x,u) 

3x(t) 

3i(x»u) 

3u(t) 

dt 

3A(tf) 

3x.(t) 

&(tf) 

3u(t) 

3A(tf) 

(2-28) 


d_ 

3A(t) 

3^(x,u, A) 

3x(t) 

3i(x.u»A) 

3A(t) 

4. 

3i(x,u,A) 

r  au(o 

dt 

,3A(tf) 

3x(t) 

,3A(tf) 

T 

3A(t) 

9£(tg), 

T 

3u(t) 

l  3A(*"£ 

(2-29) 
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c 


and: 


3G(x,u,A) 

3x(t) 

-r 

3G(x,u,j0] 

fsA(t) 

[3G(x,u,A)l[3u(t)  1_  A 

ax(t) 

3A(tf). 

3A(t)  j 

[  3u(t)  J[3A(tf)J 

where  the  final  conditions  are  specified  aa: 


3x(tf) 


and 


-  i 


(2-31) 


(2-32) 


[.an(t£) 

The  last  equation  can  be  solved  for  1 3u(t) ]  analytically  and 
substituted  into  Eq(2-28)  and  Eq(2-29)  to  eliminate  the  quantity  from 
the  equations.  The  solution  to  Eq(2-28)  and  Eq(2-29)  provides  the 
boundary  value  partials  to  be  used  to  form  the  Newton-Raphson  matrix 
in  Eq(2-21).  An  algorithm  for  Implementing  this  method  is  presented 
at  the  cad  cf  this  s actios. 

In  the  case  of  the  forward  integration  for  this  example,  the 
unspecified  states  and  adjoints  at  the  initial  point  are  iteratively 
adjusted  to  satisfy  the  specified  states  at  the  final  point.  The 


<5x(tf)  - 

which  may  be  rewritten  as: 


■ 

r  this  case  is: 

ctX<tf)j  32£<t*.) 

dw(tj)) 

3w(t0){3i(t0) 

3x(tf)j^t^ 

-1 

jx(t|) 

^ 0 ^  j 

(2-33) 


(2-34) 


where  for  tills  case: 

$x(t£)  ■  an  n-vector  of  the  deviations  of  the  final  states  from 
the  specified  values. 
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fiu/tg )  ■  a  q~ vector  of  the  necessary  changes  to  the  unconstrained 
Initial  adjolnts. 

-  an  (n-q) -vector  of  the  necessary  changes  to  the  unspec¬ 
ified  Initial  states. 

The  recursive  equation  for  this  case  la: 


“(tp) 

n+l 

“(tp) 

n 

» 

3x(tf) 

3x(tf) 

lcto> 

.  i(t0) 

T 

& 

rt 

o 

a±(t0>. 

[xf  -^(tf)] 


(2-35) 


The  perturbation  equations  for  the  terns  in  the  matrix  in  the  above 
equation  will  not  be  presented  here  but  have  the  same  fora  as  the 
equations  derived  for  the  case  of  backward  integration  presented 
earlier. 

In  order  to  satisfy  control  constraints  for  a  bounded  control 


1  4  •  *  - 

a»*4 
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cedure  is  used  (Ref.  2).  In  Valentine's  procedure,  a  new  constraint 
is  added  to  the  problem  and  is  of  the  fora: 

w(u)£.  0  (2-36) 

It  is  assumed  that  the  control  constraint  is: 


Bj,  <.  u  <  Bt  (2-37) 

or  ju  —  JB^  jO  (2-38) 

and  J9j  -  u  l  0  (2-39) 

These  two  inequalities  can  be  combined  as : 


w(u) 


(ui  -  BL1)  (B  j  -  Uj) 


(UJ  "  BJ.j)  (BTj  "  V, 


2.0 


(2-40) 
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c 


Ihis  constraint  is  added  to  the  Hamiltonian  with  a  new  Lagrange  multi- 

oiler.  ij.  T  f  t-h**  nw  Hnnrl  1  fnnl  an  ‘in  ^  hhon* 

h'  -  H  +  jyfeCu)  (2-41) 

To  satisfy  the  original  problem,  the  nev  term  must  always  be  Identic¬ 
ally  zero.  When  a  particular  control  component  is  on  one  of  the 
boundaries  then  that  component  of  w(u)  Is  zero.  If  the  control  com¬ 
ponent  Is  not  on  a  boundary,  then  the  multiplier  must  be  Identically 
zero.  A  nev  gradient  equation  1st  formed  as: 


(2-42) 


or 

H  ^  ^  wl  rn^afrA^nfv  a  nnrtfrnl 

on  a  boundary,  jj.  must  be  computed  such  that  the  gradient  equation  la 
zero.  The  fact  that: 

/w(u)  “  0  (2-44) 

Is  not  a  sufficient  condition  for: 

aw(u) 

3u 


T 

JJ  -  0  (2-45) 


G*  -  G  + 


3w(u)  1 


3u 


-  0 


(2-43) 


For  this  reason,  when  a  control  u  is  on  a  boundary, p  is  computed 

J  J 

from  the  corresponding  gradient  equation: 


-ft 

BLj  +  BTj  “  2uj 


(2-46) 
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When  the  control  Is  on  the  lower  boundary 


9H 

1*1 


Ibli 


®Tj  “  ®Lj 


and 


<  <■«  -  V  ± 


It- kn\ 


(2-48) 


which  Implies  that  Uj  Is  negative  on  the  lower  boundary.  When  the 
control  Is  on  the  upper  boundary 


btj  _  blj 


and 


YBLj 


*°  ;  <BTJ-BLJ>*° 


(2-49) 


(2-50) 


which  Implies  that  u.  la  negative.  Then  <i .  la  alwayn  nevaHve  when  the 
J  J 

control  Is  on  either  boundary.  The  term  In  the  numerator  la  positive 
because  the  original  gradient  equation,  with  no  constraints,  is  of 
the  form 

uj  -  hj(x,X)  -  0  (2-51) 

when  the  control  Is  on  the  lower  boundary: 


UJ  ■  blj 


(2-52) 


But  hj (x,X)  is  the  value  the  control  would  be  If  no  boundary  were 
present  and  is  consequently  less  than  B.,«  This  Implies  that 
Uj  -  hj (x,_X)  Is  greater  than  zero.  When  the  control  Is  on  the  upper 
boundary 


U1  mhi 


(2-53) 


But,  in  this  case,  hj(xlal)  is  larger  than  the  boundary  and  Uj-hj(x.X) 
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is  less  than  zero. 

When  the  control  is  on  a  boundary,  then  changes  in  control  at  that 
point  are  not  functions  of  changes  in  either  the  state  or  adjoint  var¬ 
iables  at  t0  or  t^  and 


auj(t) 

•  n  • 

3uj(t) 

3A(tf) 

m  u  • 

3x(tf) 

auj(t) 

3uj(t) 

3A(t0). 

“  0  ; 

3x(tfl) 

(2-54) 


at  that  point. 

When  the  control  is  not  on  a  boundary,  then' 

Uj  -  0  (2-55) 

and  the  control  is  computed  as  If  the  boundary  on  control  did  not 
exist. 

Au  algorithm  Lot  implementing  the  second  variation  procedure  with 
control  inequality  constraints  and  using  backward  integration  follows. 
The  method  for  forward  integration  is  similar. 

Step  1.  An  initial  set  of  estimates  are  made  for  the  ter¬ 
minal  adjoint,  variables. 

Step  2.  The  control  is  computed  at  tf  from  the  gradient 
Eq(2-43).  If  it  la  within  bounds  go  to  Step  3;  otherwise,  go  to  Step 
5. 


Step  3.  The  state  and  adjoint  equations -are  integrated 
backward  for  one  increment  of  time. 

Step  4.  The  control  is  computed  at  this  point  from 
Eq(2-43).  If  it  is  within  bounds  go  to  Step  6;  otherwise  go  to 
Step  5. 
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Step  5.  Let  the  control  equal  BT  or  and  compute  from 
Bq(2-47).  Set  [8uj (t)/31(tf) ]  equal  to  tero  and  go  to  Step  7. 

Step  6.  Set  p  -  0  and  continue. 

Step  7.  Solve  the  perturbation  equations  backward  In  time 
to  the  same  point  as  the  states  and  adjoin ts. 

Step  8.  If  the  initial  time  has  been  reached  go  to  Step  9; 
otherwise  go  to  Step  3. 

Step  9.  Compute  new  estimates  of  the  terminal  adjoints 
from  Eq(2-22)  by  inverting  the  nxn  matrix  and  go  back  to  Step  2  and 
continue  until  the  boundary  conditions  converge  to  within  a  tolerance 
of  the  specified  values. 

Suboptimal  Closed-Loop  Control  Methods 

The.  basic  purpose  of  using  a  closed- loop  control  in  conjunction 
with  an  optimal  problem  is  to  form  a  practical,  near  optimal,  guidance 
law  for  use  in  a  real  time  adaptive  controller.  If  this  suboptimal 
law  is  used  with  perfect  measurements ,. the  feedback  elements  derived 
herein  will  satisfy  the  final  boundary  conditions  established  in  the 
optimal  solution.  The  solution  of  the  open-loop  optimal  control 
problem  determines  a  nominal  trajectory  and  control  set  which  satisfy 
the  boundary  conditions  at  both  ends  of  the  trajectory.  However, 
when  the  optimal  control  is  applied  as  a  forcing  function  during  a 
mission,  the  actual  initial  conditions  may  not  be  identical  to  the 
nominal  initial  conditions  for  the  optimal  trajectory.  Also,  due  to 
perturbations  along  the  trajectory,  the  actual  trajectory  may  deviate 
from  the  nominal  optimal  trajectory  by  a  small  amount.  For  a  mission 
such  as  entry  into  the  earth's  atmosphere,  there  still  exists  a 
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requirement  to  reach  a  particular  target  with  a  set  of  specified  final 
states.  An  optimal  closed- loop  control  scheme  corrects  for  these  dev¬ 
iations  by  continually  computing  a  new  optimal  trajectory  to  arrive  at 
the  specified  boundary  conditions.  However,  the  computation  of  an 
optimal  control  at  every  point  in  time  to  satisfy  the  final  boundary 
conditions  requires  the  complete  solution  of  the  original  optimal 
control  problem  from  the  actual  initial  conditions  at  each  point  to  the 
final  point. 

An  alternative  to  solving  the  entire  problem  at  each  point  is  to 
define  a  linear  functional  relationship  between  small  deviations  in 
the  states  from  the  nominal  and  small  changes  in  control  from  the 
nominal  to  achieve  the  desired  final  boundary  conditions.  This  func¬ 
tional  relationship  is  only  valid  in  a  linear  range  about  the  nominal, 
but  allows  a  practical  comoensAtlon  for  small  errors.  This  functional 
relationship  is  given  in  reference  1  as: 


(2-56) 


where  dju(t)  is  the  required  change  in  control  vector  to  correct  for 
small  errors  <5x(t)  in  the  state  vector.  The  matrix  is  computed  from 
quantities  available  from  the  original  solution  of  the  optimal  open- 
loop  control  problem.  It  should  be  noted  that  the  change  in  control 
generated  does  not  immediately  correct  the  trajectory  to  the  nominal. 
Instead,  it  allows  the  states  to  follow  a  neighboring  path  to  the 
specified  final  conditions.  This  is  illustrated  in  Figure  1.  In  order 
to  accomplish  this,  the  matrix  [3u(t)/  3x(t)]  is  computed  backward  in 
time,  relating  all  variables  to  the  final  specified  conditions.  This 
matrix  is  referred  to  as  the  feedback  gain  matrix. 
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Figure  1.  Neighboring  Trajectory  for  Linear 
Feedback  Scheae 

The  block  diagram  In  Figure  2  llluatratea  an  axampla  of  the 


Implementation  of  a  linear  feedback  echeme. 


yfd-)  stored 
in  computer 


Physical 
System  Dynamics 


% 


Figure  2.  Block  Diagram  for  Implementing  Feedback  Galna 
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The  quantity  which  must  b«  derived  is  the  rslstionshlp  for  the 
feedback  gain  matrix.  The  equation  for  defining  this  matrix  as  derived 
by  Johnson  in  Reference  1  la: 

au(t)  _  _  ac(t)  -1  3G(t)  +  3G(t)  ax(t)  a«(t)  ,, ^ 

ax<t)j "  *"[  au(t)J  |^ax(t)j  ax(t)  ax(tf)  ax(tf) 

where  G(t)  la  the  gradient  and 


3X(t)  3x(t) 

ax(t,)  *nd  ax(t2) 


are  computed  by  solving  the  perturbation  equations,  Eq(2-28),  Eq(2-29), 
and  Eq(2-30)  in  the  previous  section.  This  equation  la  only  valid  for 
the  case  in  which  the  states  are  specified  at  both  ends  of  the  trajec¬ 
tory.  It  is  not  practical  to  force  the  change  in  control  to  be  a 
function  of  errors  In  the  nrHn-J cannot 
be  generally  computed  during  an  actual  entry  mission. 

The  algorithm  for  computing  the  feedback  gains  simply  consists  of 
integrating  the  three  perturbation  equations  backward  in  time,  substi¬ 
tuting  the  quantities  into  Eq(2-58),  and  storing  the  resulting  gain 
values  at  aach  point  in  tima.  It  should  be  noted  that  these  gain 
values  cannot  be  computed  from  forward  integration  because  the  final 
specified  boundary  conditions  do  not  directly  enter  into  the  lntegra- 


Qptlmal  State  Estimation 

The  optimal  state  estimator  used  in  this  thesis  is  the  discrete 
Kalman  filter.  The  Kalman  filter  is  a  minimum  variance,  unbiased, 
linear  estimator.  The  classical  least  squares  estimator  is  a  special 
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case  of  the  Kalman  filter.  Whan  no  plant  nolsa  la  praaant  and  tha 
plant  differential  equationa  are  uncoupled,  the  Kalman  filter  and  the 
unweighted  least  aquarea  estimator  are  equivalent.  The  discrete  Kal¬ 
man  filter  is  expressed  ss  s  set  of  recursive  equations  just  as  the 
least  squares  estimator  can  be.  A  derivation  of  the  filter  equations 
will  not  be  praeented  here  but  la  well  described  In  Kef.  3,  12, 
and  13. 


When  the  plant  and  observation  noise  are  independent,  white 
Gaussian  vectors  and  the  plant  and  observation  equations  are  linear, 
the  Kalman  filter  provides  tha  optimal  Bayesian  estimate  of  the  state 
vector.  When  the  noise  is  not  Gau  elan,  the  Kalman  filter  still 
provides  the  optimal  linear  estimate  of  tha  state  vector.  When  the 
plant  or  observation  vector  equationa  are  non-linear,  the  optimal 


^  ^  d.  *•» «»  •  «  a  J  h,  • 
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lved  for  each  problem  individually  whether  the  noise  vectors  are  white 
Gaussian  vectors  or  not.  For  most  non-linear  problems,  the  most 
practical  approach  is  to  linearise  the  equations  and  apply  the  Kalman 
filter  to  the  small  deviations  within  a  linear  range.  The  following 
basic  procedure  is  generally  followed  in  applying  the  Kalman  filter  to 
a  non-linear  set  of  equationa: 

A.  The  non-lineA--  plant  and  observation  equations  are 
linearized  either  about  a  precomputed  noulnal  set  of  values  of  the 
state  vector  or  about  the  current  optimal  linear  estimates  of  the  state 
vector.  In  this  thesis,  only  linearization  about  the  current  optimal 
linear  estimates  are  considered. 

B.  The  Kalman  filter  is  applied  to  the  linear  set  of 
deviations  from  the  nominal  or  estimated  states  to  form  the  optimal 
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estimates  of  the  deviations » 

C.  The  optimal  estimates  of  the  deviations  are  added  to 
the  precomputed  nominal  states  or  to  the  best  estimate  of  the  states 
given  nil  information  except  the  current  observation  vector  values • 

Thie  provides  the  optimal  linear  estimate  of  the  state  values. 

The  system  model  (l.e.,  plant  model,  message  model)  equation  and 
the  observation  equation  are  defined  a a  follows: 

x  “  i(x,I,u>  (2-58) 

i  -  j &(x,T])  (2-59) 

whera  x  •  stats  vector  to  be  estimated 

jr  «  observation  vector 
u  -  control  vector 
£  ■  Gaussian  white  system  noise 

n  ■  nnna«4  on  4  fa  HCf.CC 

In  order  to  establish  e  discrete,  recursive  set  of  equations,  the 
system  and  observation  equations  are  written  in  the  following  discrete 
form: 

~inH~  "  "  I(* i -ii , Jit)  (2-60) 

or  *i+l"  2£i  +  H*I(xi»ii»iii)  4  Z(it»ii»iii)  (2-61) 

end  JLf  “  4£(25i»Hi)  (2-62) 

where  U  is  the  integration  step  size  for  a  linear  approximation.  When 
the  discrete  system  and  observation  equations  are  expanded  in  a  Taylor 
aerlas  about  the  optimal  estimates  1^*4^  and  Hi  the  following 
equations  are  obtained: 
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li  "  G(*i>Q) 


(*i  ~  *i>  + 


Xj»Xj 


aGCxj.Ju) 


am 


m 


(2-69! 


TW-0 


The  Kalman  filter  may  be  used  with  theae  equations  In  two  ways. 
Both  are  equivalent.  One  method  Involves  applying  the  linear  filter 
equations  to  the  deviations  from  the  optimal  estimate  to  find  the  op¬ 
timal  estimate  of  the  deviations  and  then  adding  the  result  to  the 
best  estimate  of  the  states.  The  other  method  consists  of  modifying 
the  filter  equations  to  include  the  best  estimates  of  the  states.  The 
latter  procedure  will  be  followed  here.  The  best  estimate  of 
given  only  2.1  and  all  preceding  observations  is  denoted  x(i+l| 1) .  The 
best  estimate  may  be  computed  by  solving  Eq(2-58)  as: 

x<  j  1 1)  -  f  [x(  i  I « .0 .  uj  ]  ilj  <1+1  <2-70: 

where,  the  integration  begins  at: 

x(l|l)-ii  <2-7i: 

and  ends  when  j  -  1+1. 

The  recursive  relations  to  be  computed  for  the  filter  are  as 
follows : 

A.  Kalman  Gain: 

Ki+l  -  T (1+1 1 1) Gx<  1+1) T [ Gx(  1+1) r ( 1+1 1 1) Gx(  1+1) T  +  w^l'1  (2-72: 

B.  Best  Estimate  Of  Error  Convarlance  Given  an: 

r(i±i|«  -  f^f*  +»t  (2-7* 

This  quantity  may  be  computed  by  propagating  over  small  Intervals  with 
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the  state  transition  matrices  Fx  and  F  at  the  same  tine  x(.<+l|l) 
is  computed. 


r.  Wrrnr  fnvnWiini'#t  Clvim  Vint 

- "  ZZ - ~~  - 

ri+l  -  U  -  Ki+iGx(i+l)  r(i+l|i)tl  -  Ki+i^d+l)]1 

+  K1+1W1+1K1+1T 


(2-74) 


D.  Optimal  Estimate: 


*i+l  “  *<1+1!  V  +  Ki+1  *1+1  "  G[x(i+l|  i)  ,0] 


(2-75) 


In  the  four  sets  of  equations  above t 

^i+1  ”  Kalman  gains  at  measurement  point  (i+1). 

'  r(i+l|i)  -  Best  estimate  of  error  covariance  at  measurement 

point  (1+1)  given  Xi* 

•  ( 

ri+1  ■  Error  covariance  at  measurement  point  (1+1) 

given  jq+1 

x(i+l|i)  -  Best  estimate  of  state  vector  at  measurement  point 
(i+1)  given  y^. 

-  Optimal  estimate  of  state  vector  at  measurement 
point  (i+1)  given  Xi+1* 


a*i+l 


*1+1  "  *(i+i|i) 


a*i 


*1  ■  *i 


v1  -  ECvjVi1) 


where 


31(xi»Ci*Jii> 
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■  ■:  ■  ■-  -*  ,  >V ,  ■ 

r  •  f  - 


y  *  v . 


c 


wi "  E<Hi«i> 
where 

2.1 


3_H. 

W  1  J 


%"0 


*i 


A  specific  algorithm  for  implement lag  the  Kalman  filter  will  be  pre¬ 
sented  in  Chapter  VI  on  Optimal  Implicit  Guidance.  These  equations 
will  form  the  basis  for  the  optimal  estimation  in  the  entry  guidance 
law  formulation. 


Initial  Estimate  Of  The  Error  Covariance  Matrix.  In  this  sec¬ 
tion!  a  set  of  equations  will  be  derived  for  computing  the  initial 
error  covariance  matrix,  r i*  t  when  certain  assumptions  are  satisfied. 

It  will  be  shown  that  it  is  not,  in  general,  necessary  to  guess  an 
initial  estimate  To*  It  stay  be  computed  Instead.  It  is  assumed  here 
that  the  problem  begins  at  the  first  measurement  point.  Theoretically, 
for  a  linear  system  and  linear  observations,  given' a  long  enough 
sequence  of  measurements,  any  estimate  of  the  initial  error  covariance 
matrix  should  produce  the  same  steady  state  errors.  In  actual  practice 
with  non-linear  systems,  given  non-linear  observations  and  a  finite 
sample  time,  this  does  not  occur.  This  is  illustrated  by  Elliott  and 
Fillatreau  in  Ref.  17.  They  investigated  the  effects  of  variations  in 
the  initial  estimate  of  the  error  covariance  matrix  in  conjunction 
with  an  orbit  determination  problem. 

The  steady  state  errors  are  sometimes  critically  dependent  on  the 
Initial  estimate  of rQ,  If  the  estimate  is  much  smaller  than  the  in¬ 
itial  observation  -loise  matrix,  W(0) ,  the  variables  may  not 
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reach  steady  state  values  at  all.  If  the  estimate  is  larger  thanW(O), 
the  variables  in  soma  cases  reach  steady  state  values  in  a  much  shorter 
time  than  for  lover  estimates.  When  the  estimate  la  much  larger  than 
W(Q) ,  the  reference  cited  above  indicates  a  tendency  for  the  variables 
also  not  to  reach  steady  state  values  as  in  the  case  in  which  the 
estimate  is  much  smaller  thanW(O). 

The  equation  for  the  optimal  estimate,  Eq(2-75),  is: 

ki+l  -  I U  +  Ki+1 &1+1  -  G(x(i+1 | i) ,0) )  (2-76) 

For  the  initial  observation  point  this  equation  becomes! 

*0  “  i°  +  K0  b*  ~  G(x°,0)]  (2-77) 

where  x°  is  defined  as  tha  bast  estimate  of  the  initial  state  vector 
given  no  observations,  and  x^  is  the  optimal  estimate  of  the  initial 
state  vector  given  the  initial  observation.  Thus,  can  be  based 
nnly  on  s  priori  information.  The  following  assumptions  arc  aov  made; 

A,  x°  is  an  arbitrary  a  priori  estimate  of  the  initial 
state  vector. 

B.  KQ  is  an  arbitrary  a  priori  estimate  of  the  initial 
Kalman  gain  matrix.  The  value  of  KQ  will  be  determined  by  the  con¬ 
fidence  level  of  the  a  priori  x°  as  compared  to  the  confidence  level 
of  the  observation  vector  as  given  by  W(o) . 

The  equation  for  the  initial  Kalman  gain  from  Eq(2-72)  is: 

K0  -  r°Gx(0)T(Gx(0)r°Gx(0>'r  +  W(O))-1  (2-78) 

where  r°  is  defined  as  the  best  estimate  of  the  initial  error  covari¬ 
ance  matrix  given  no  meaaurments.  The  optimal  estimate  of  the  initial 
error  covariance  matrix  given  the  first  measurement  will  be  denoted  r°. 
Solving  Eq(2-78)  for  gives: 

r°  -  [I  -  K0Gx(O)]"1K0W(O)Gx(O)-* 
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The  equation  for  the  optimal  estimate  of  the  Initial  error  covariance 
matrix,  r0 ,  from  Eq(2-74)  le: 

r°  -  [I  -  K0Gx(0)]r°[I  -  KoG^O)]1  +  K^WCOKo1  (2-Wi> 

Substituting  Eq(2-79) : 

r0  ■  [i  -  k0gx(o)][i-k0gx(o)]"1k0w(o)gx(o)-t[i-k0gx(o)]t+k0w(o>kJ(2-81) 


or! 


r0  -  K0W(O)Gx(O)“T[I  -  K0Gx(0)]T  +  K0W(O)  Kj  (2-82) 

Simplifying  this  equation  gives: 

ro  -  K  W(0){Gx(0)"T  -  Gx(0)"T[K0Gx(0)]T  +  (2-83) 

or: 


rQ  -  KflW(0)Gx(0)-T 


(2-84) 


where: 


*  *  -  •#  —  — —  -  — 


Gx«» 


th 


—  —  g  m  «mwb 


(2-85) 
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matrix,  rQ,  once  Xq  has  been  computed  using  arbitrarily  selected 
values  for  x°  and  Xq. 

A  special  solution  of  this  problem  will  now  be  considered  In 
which: 

K0  4  1  (2-87) 


and: 


x°  -  G(x°  ,0)  (2-88) 
In  the  case  in  which  the  observation  vector  Is  a  linear  function  of  the 
state  vector  or: 


■  *i  +  Hi  (2-89) 

.  Eq(2-88)  is  satisfied  by  any  value  of  jx°  and  Eq(2-85)  and  Eq(2-84)  give 

Gx(0)  -  I  (2-90) 
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and: 

r0  -  W(0)  (2-9 X) 

Ilia  anaa^al  ana  a  df-Cfi^.d  «KavA  4  M  f)lA  AAA  uK^h  u{1  1  Ka  AnnHfl^  ill 

tha  optimal  implicit  guidance  aimulation  for  this  theaia. 

Eq(2-84)  provides  a  relationship  which  allows  the  following 
Initial  estimate  algorithm  to  be  constructed: 

Step  1.  By  some  arbitrary  method,  select  an  xa  and  Kq 

such  that: 


^  +K0CXo  -  G(*  ,0>]  (2-92) 

provides  the  desired  relationship  for  x0  (such  as  *PmQ  and  Kg "I) 


Step  2.  Evaluate : 


0,(0) 


G(x,0) 


T 

JMf 


(2-93) 


and  invert  the  matrix  to  form  G^0)“T. 

Step  3.  Compute  the  initial  error  covariance  matrix 


from: 

r0  -  K^OG^O)"1  (2-94) 

This  method  of  estimating  the  initial  error  covariance  matrix 
and  a  further  development  of  the  case  to  which  it  la  applied  are 
presented  in  Chapter  VI.  The  following  chapter  shall  develop  the 
model  equations  for  the  optimal  open-loop  control  problem. 


B 


« 


L 
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III.  Derivation  Of  Basic  Entry  Equations 

In  this  chapter  the  basic  equations  of  motion  for  an  entry  body 
will  be  established  along  with  the  equations  for  the  position  sensi¬ 
tivity  coefficients  to  be  minimized.  The  lumped  parameter  uuuti'ol  con 
taining  parameters  for  vehicle  configuration  and  attitude  will  be 
defined  and  the  characteristics  and  effects  of  the  atmosphere  will  be 
discussed. 

State  Vector  Equations  Of  Motion 

Figure  3  below  shows  the  geometry  and  sign  conventions  to  be  used 
for  the  entry  problem.  The  motion  of  the  entry  vehicle  can  be  descri¬ 
bed  by  a  set  of  four  coupled  differential  equations.  The  solution  of 
these  equations  provides  position  with  respect  to  the  center  of  the 
earth  (r^,  velocity  (v).  areal  velocity  (w)  ,  and  range  angle  (o) .  Fnr 


Figure  3.  Entry  Geometry 
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this  problem  a  transformation  of  th*  Independent  variable  la  mada  such 
that  the  range  angle,  c  ,  la  the  Independent  variable  Instead  of  time, 
this  reduces  the  nunber  of  equations  to  three.  A  transformation  la 
also  made  from  the  areal  velocity,  v,  to  flight  path  angle, Y. 

The  equation*  of  morion  for  a  planar  trajectory  with  reapact  to  a 
spherical,  non-rotating  earth  with  ablation  effects  neglected  aa 
developed  by  Johnson  In  reference  1  are: 


i-  -X 
r 


V  •  -  -$p3-D 


w  ■  v2-  ^  -Dtcos(y)  +Lrsin(Y) 
r 

A  m  va*n(Y) 


(3-1) 

(3-2) 

(3-3) 

(3-4) 


Where  r  ■  displacement  with  respect  to  the  center  of  the  earth, 
v  ■  velocity. 

Y  ■  angle  of  the  velocity  vector  measured  from  the  local  verti¬ 
cal.  This  is  always  a  positive  quantity, 
o  -  range  angle.  The  range  angle  la  a  positive  quantity  and 
assumed  to  be  a  monotonic  function  of  time, 
w  -  areal  velocity.  This  quantity  Is  defined  by  the  relation 
w  ■  rvcos (y) • 

D  “  drag.  The  drag  is  dependent  on  the  coefficient  of  drag  and 
the  square  of  the  velocity  and  is  normalized  to  be  drag 
force  per  unit  mass. 

L  -  lift.  The  lift  is  dependent  on  the  coefficient  of  lift  and 
the  square  of  the  velocity  and  is  normalized  to  be  lift 
force  per  unit  mass. 
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p  ■  universal  gravitational  constant  times  the  mass  of  the  earth. 

This  quantity  is  the  eerth's  gravitational  constant. 

The  areal  velocity,  v,  is  used  to  improve  the  stability  of  the 
scdsl (  for  Isrjs  viriiticr.s  iti  y -  Since  only  ■  thirty  degree  wu> vi twin 
variation  in  yis  expected,  it  is  not  necessary  to  use  the  areal  veloc¬ 
ity.  Also  the  computation  of  w  onboard  an  entry  vehicle  would  require 
a  continual  computation  of  (cosy  )  throughout  the  flight.  The  quantity 
y  ,  is  available  onboard  from  the  glmbal  resolvers  in  the  Inertial 
Measuring  Unit  (IMU).  For  these  reasons  a  new  set  of  equations  will  be 
derived  to  rnplace  w  with  y  . 

The  relationship  between  w  and  y la : 

w  ■  rvcosy  (3-5) 

When  the  time  derivative  of  w  is  computed,  the  result  is: 

$  ■  fcveosy  +  rvcosy  -  rvfsiny  (3-6) 

Substituting  the  equations  for  b  and  b  from  Eq(3-1)  and  Eq(3-2)  gives: 

*  -  wvcoaY  -  HggPP.Y,  -  rDcosy  -  rvVsiny.  (3-7) 

r  vrz 

or 

Veiny  -  -  ^2S1  _  *-  (3-8) 

r  *  v*r3  v  rv 

Substituting  for  *  from  Eq(3-3)  and  for  w  from  Eq(3-5),  combining 

terms  and  dividing  by  slny  gives: 

•  „  _  veiny  _  usinv  _  L 

r  rzv  v 

A  new  set  of  state  equations  may  now  be  written  as  follows: 

b  ■  vcosy 

*  -  -  -HSSS1  -d 

xt 


(3-9) 

(3-10) 

(3-11) 
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•  _  veiny  +  ualnV  _  £. 
T  r  r*v  v 


•  u«4fiv 


(3-12) 


(3-13) 


Th«  independent  variable  In  the  above  equatlona  la  time,  which  la 
monotonic  and  naaaurabla.  It  la  aaauaed  that  ranga  angle  la  alao  mono- 
tonic.  To  alnplify  the  modal  equatlona,  the  Independent  variable, 
time,  la  replaced  with  ranga  angle  cr.  It  la  important,  however,  to 
compute  time  versus  range  angle,  to  be  able  to  revert  to  a  time  baae 
at  tha  conclusion  of  the  problem. 

Dividing  Eq(3-ll)  through  Eq(3-13)  by  ft  yields  the  following  set 
of  dynamic  equatlona: 


r 

tany 


(3-14) 


v'  A  ibi  m  -  U-  _  rp 
*  da  rvtany  vsiny 


(3-15) 


r'4^-- 

da 


1  + 


_JU _ tk. 


rv‘ 


fZ8iny 


(3-16) 


To  determine  the  time,  t,  at  any  point  on  the  trajectory,  the  recip¬ 
rocal  of  Eq(3-13)  can  be  used,  which  is: 

l'  4  H  ■  v.In(YT  °-1,) 

Introducing  state  space  notation  to  the  model  equations  derived 
above  gives: 

xT-  [xitx2,x3]  ^  (r,v,y]  (3-18) 


/'.nr/vvJin-tn 
•**•»  »  —  -- 
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_d 

8x,<Of) 

_ 

3x,(Of)  8x, (Of) 

3x,(0f) 

do 

8x(c) 

do 

axjCo)  *  8x2 (c) 

’  8x3(0) 

Based  on  this  and  Eq(3-23): 


[  x4,  x5,  x6]  -  -  [  x4,  x5,  x6] 


r  3fi  3fi  3fii 


8xi 

8x2 

Hv« 

8f2 

CM 

3xi 

8x2 

8x3 

®f  3 

*1 

n— 

8x3. 

(3-33) 


(3-34) 


where  the  partlals  are  defined  In  Eq(3-24)  through  Eq(3-32). 


Lumped  Parameter  Control  Equations 

The  design  of  a  passive  control  restricts  the  control  available 
to  atmospheric  forces.  In  this  section  the  equations  for  drag  and 
•lift  shall  be  functionally  separated  into  two  parts.  The  first  part 
shall  contain  parameters  related  to  vehicle  configuration  and  attitude 
and  shall  be  defined  as  a  lumped  parameter  control.  The  second  part 
shall  contain  parameters  related  to  mach  number. 

The  equations  for  lift  and  drag  are: 


u  _  AC6(“)  r  pv2cpa 

m  2 

(3-35) 

L  _  SCL<a>  ,  Pv2cos 
m  2 

(3-36) 

where  p  -  density 
v  *  velocity 

CpS  -  coefficient  of  pressure  at  the  stagnation  point,  a  func¬ 
tion  of  mach  number 

A  ■■  cross-section  area  perpendicular  to  velocity  vector 
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S  •  plan fora  area  perpendicular  to  both  tha  velocity  vaetor 
and  tha  radlua  vaetor 

a  *  aaaa  (a  constant) 

Cl (a)  *  a  draff  coafficlant  wh(rh  (a  a  function  of  an«1a  nf 

attack,  a 

C^(a)  ■  a  lift  coefficient  which  la  a  function  of  angle  of 
attack,  a 

Two  controls  will  be  defined  ss  follows: 


“dA. 


ACfc(a) 


(3-37) 


(3-38) 


The  values  which  the  controls  are  allowed  to  take  on  are  limited 
by  practical  considerations  of  vehicle  design,  thus  the  lumped  control 
problem  is  a  bounded  control  problem.  The  boundaries  are  as  follow . : 
2.76  x  10” 5  <  <  2.22  x  10“3 


-1.6  Ujj  £  uL  £  1.6  uD 

These  boundaries  on  the  controls  correspond  to: 


250  £  _^£_  £  20,000  Kg/meters2 
CqA 


-1.6  £  L/D  £  1.6 


Substituting  the  control  relations  Into  Eq(3-35)  and  Eq(3-36)  gives: 


pvzC 


r“  ■« 


(3-39) 


PVZC, 


22. 


2  -»L  (3"40> 
The  term  CpS,  the  coefficient  of  stagnation  pressure  In  the  above  equa¬ 
tions  is  a  function  of  mach  number,  which  in  turn  is  a  function  of 
altitude,  pressure,  density,  and  velocity.  This  term  is  derived  using 


■\2 
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the  following  assumptions: 

A.  Perfect  Caa 

B.  Isen tropic  Flow 

C.  Adiabatic  Flow 

X).  'the  shock  ia  locally  normal 

In  Fig.  4  below,  P0l  and  P02  are  the  stagnation  pressures  at 
points  1  and  2;  P1  and  P2  are  the  pressures  at  the  two  points,  and  M 
ia  mach  number.  Using  assumptions  A,  C,  and  D,  the  following 
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relations  *rc  obtained  (Ref.  15) s 

Cr 


:P»  "  S02  "  ^  ^  _1l 


‘  l  * 


and 


Id 

POI 


Y+l’ 


1 

'  y^t 

(Y+1)H? 

(Y-DM2+2 

A 


(3-41) 


(3-42) 


From  assumption  B: 


*21  . 

e, 


A 


(3-43) 


Multiplying  Eq(3-42)  and  Eq(3-43)  together  and  substituting  the 
result  for  the  quantity  P02/Pj  in  Eq(3-41) 


ps  Y 


Y+l  . 

1 

7=T 

Y+l 

2y 

2 

#r_  2 

TSf 


(3-44) 


For  air,  Y-1.4,  which  when  substituted  into  Eq(3-44)  provides 


the  final  relation  for  C  . 

ps 


238.87872M5 

1 

7(2.8M2~0.4)2 

2(2.8Mz-0.4) 

»  • 

t 


(3-45) 


The  above  equation  is  valid  for  velocities  above  Mach  2.  The 
density  p  is  Interpolated  numerically  from  the  1959  ARDC  Standard 
Atmospheric  Tables. 

Substituting  the  expressions  of  Eq(3-39)  and  Eq(3-40),  for  drag 
and  lift  respectively,  into  the  state  equations,  Eq(3-19)  through 
(3-21)  the  final  state  equations  are: 


V* 

x\ 


__xO_ 

tanx- 


k  Fi 


(3-46) 
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/  *  L  *t\ 
\  **  -t*  f 

(3-48) 

(3-4S) 

(3-50) 

(3-51) 


These  are  the  basic  model  equations  to  be  used  In  the  optimal 
control  problem  to  be  discussed  In  the  following  chapter. 


45 


GGC/11/ 70-10 


IV.  The  Maximum  Principle  Applied  To  The  Entry  Problem 

In  Chapter  II,  it  vu  established  that  an  arbitrary  critarion 
could  ba  satlafiad  along  with  a  aat  of  non-llnaar  plant  aquations  pro¬ 
vided  a  sat  of  nacassary  optimal  conditions  ara  mat.  In  Chapter  III, 
mathematical  modal  equations  ware  derived  for  the  entry  problem.  The 
quantities  to  ba  used  in  the  criterion  function,  the  position  sensi¬ 
tivity  coefficients,  were  also  explained  in  concept  and  derived  in 
Chapter  III.  In  this  chapter,  the  criterion  function  is  established 
for  the  entry  model  along  with  a  set  of  necessary  optimal  conditions 
such  that  the  model  equations  ara  satisfied. 

Criterion  Function 

The  criterion  function  is  a  scalar  valued  function  which  contains 
all  the  elements  to  be  minimized.  For  this  problem,  the  criterion 
function  must  contain  a  convex  function  of  the  errors  at  the  terminal 
point  due  to  errors  at  any  point  of  the  entry  trajectory.  The  posi¬ 
tion  influence  coefficients,  as  developed  in  Chapter  III,  describe 
these  errors,  and  therefore  serve  as  elements  of  the  criterion  func¬ 
tion.  The  remaining  elements  of  the  function  are  the  lift  and  drag 
control  quantities.  Since  control  requires  an  expenditure  of  energy, 
it  is  generally  minimised  with  respect  to  zero.  In  this  problem,  the 
control  vector  is  in  a  lumped-parameter  form  containing  parameters 
defining  the  vehicle  configuration  and  a  variable  function  of  angle  of 
attack,  It  is  desirable  to  minimize  the  control  quantities  about 
nominal  mean  values  to  satisfy  practical  limitations  on  the  vehicle 
configuration.  This  allows  for  control  deviations  from  the  nominal 
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design  value*  a*  well  a*  determine*  the  optimal  control  baaed  on  tra¬ 
jectory  consideration*.  The  criterion  function  is: 


j  -  i  ff  [  y7RY  +  VTSV  ]  do 
2  o.  ~  ~  ~  “ 


(4-1) 


where: 


U  -  M  ' 
D  uD 

V 


and  R  and  S  are  positive  definite  diagonal  weighting  matrices. 


(4-2) 


Necessary  Conditions  For  Optimality 

The  necessary  conditions  for  optimality  are  established  by  form¬ 
ulating  the  Hamiltonian  and  applying  the  appropriate  conditions.  The 
State  vector _X  includes  the  «t*te  of  the  entry  mn^nl  ond  rho  pnairinn 
Influence  coefficients  derived  in  Chapter  III.  The  minimization 
problem  la  stated  as: 


minimize  J 


1  / 

2  ;o 


Of 


$dd  subject  to  x'  ■  „f(x,u) 


(4-3) 


The  state  equations  from  Chapter  III  are: 


x 


X 


X 


I 

1 


I 


2 


I 


3 


M 


tan  x. 


_ m _  _  Fail,  u 

x^XjtanXj  x^sin  x3  D 


-1  +  3  “  -  2-  - UT 

X1X2  X2  Bil*  X3  ^ 


L  F 


1 


4  F 


2 


(4-4) 

(4-5) 


(4-6) 
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<  “ 

-x  3Fl  - 

4  ^7 

x  3F2  - 

Xs .*7 

xe 

3?3 

IXj 

L  f4 

(4-7) 

3F1 

3F, 

3F, 

„  1  , 

A  V 

//._G\ 

--5 

"4  3x2 

"5  3x2 

"6 

3x2 

-  5 

\-» 

3Fj 

"X4  ixj  " 

3F, 

3x7 

(4-9) 

X*  = 
x6 

xs  ^7 " 

x6 

A  *6 

The  complete  equations  of  the  necessary  conditions  for  optimality 
are  not  developed  here,  but  are  included  ip  Appendix  A.  For  the  dev¬ 
elopment  in  this  section,  a  set  of  equations  only  at  a  functional 
level  is  given. 

Given  the  formulation  of  the  minimization  problem  in  Eq(4-3) ,  the 
Hamiltonian  is  as  follows: 

£ 

H  ■  1  +  i:  >.rF|  (A-1™ 

i-1  1 

where  «  state  derivatives 

A^  ■  adjoint  variables 

With  the  Hamiltonian  defined,  the  four  sets  of  necessary  condi¬ 
tions  may  be  stated  for  the  entry  problem  by  referring  to  Eq(2~10) 
through  Eq(2-13)  in  Chapter  II. 

The  State  Squat Ions.  The  yte  equations  have  been  stated  in 
Eq(4-4)  through  Eq(4-9) ,  and  are  given  by: 


3H 

3Ai 


(4-11) 
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Table  II,  Boundary  Conditions  For  The  Entry  Problem 


a0 

Of 

xi 

(i  -  1.2,3) 

specified 

specified 

\ 

"I 

(i  -  1,2,3) 

unconstrained 

unconstrained 

XJ 

(j  -  4,5,6) 

unconstrained 

specified 

(1,0,0) 

(j  -  4,5,6) 

0 

unconstrained 

When  all  four  conditions  are  satisfied,  the  optimal  problem  Is 
solved.  To  accomplish  this,  the  numerical  procedures  described  In 
Chapter  II  are  us^d.  The  application  of  the  algorithms  to  the  optimal 
entry  problem  Is  discussed  in  the  following  section. 

Second  Variation  Equations 

The  complete  perturbation  equations  for  the  second  variation  pro- 
cedure  are  not  presented  here  but  are  Included  in  Appendix  A.  For  the 
entry  problem,  It  is  assumed  that  position,  velocity,  and  flight  angle 
are  specified  at  the  initial  point  and  the  sensitivity  coefficients 
are  unspecified.  It  is  also  assumed  that  all  six  states  are  specified 
at  the  terminal  point.  Thus  for  backward  Integration,  iterative  esti¬ 
mates  are  made  on  all  six  terminal  adjoint  values.  The  Newton-Rephson 
equation  for  this  case  is: 
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P_.  /_  si  [3  xi  (o0)  8xt  (0O)_  3xi  (on)  3xi  (ap)  3xt  (op)  3xi  (oq) 

0  l'  aXiCap  3X2(a£)  3Xj(af)  3X4(of)  3X5(af)  3X6(of) 

^  (n  ^  3x2(an)  3x2 (op)  3x2(a0)  3x2(g0)  3x2(oq)  3x2(ao) 

axpop  3X2(of)  3Xg(of)  Ms{0f)  »A6(o£i  | 

3v-i(riQ)  3xo(op)  3xg(jn)  8x3(03)  3x2(gp)  3x3(g n)  I 
liX3''  3X^(Or)  3X2(o£)  3X3(o£)  3X4(o£)  3X3(a£)  3X6(a£) 


3X4(00)  3X4(ao)  3X4(oq) 


3X1 (af)  3X2(of)  3X3(of)  3X4(o£)  3X5(o£)  3X6(o£) 

3X5(03)  3Xj(0q)  3Xs(Oq)  SXgCOp)  3X5(Oq)  3X4(Oq) 


6Xc(a<r) 


3Xj (oo)  3X2(oo)  3X3(03) | 8x4(00)  8x5(00)  3x6(°o) 

3x2(Qf)  3x2 (of)  3x2 (0f) | 3x2 (of)  3x2 (gf)  3x2 (of) 
8X^03)  3X2(o0)  3X3(oo)|3x4(0o)  3x5(oq)  3xg(o0> 

3x3(gf)  3x3(0f)  3xq(gf) | 3x3(gf)  3x3(gf)  3x3(gf) 
SX^Oq)  3X2(o0)  3X3(o0)  |3x4(a0)  3x5(o0)  3x6(o0) 

3x4(0f)  3x4(gf)  3x4(of) I 3x4 (gf)  3x4(gf)  3x4(gf) 

3Xj(Oq)  3X2(Oq)  3X3(0q) I 3x4(0q)  3Xg(0g)  3Xg(0Q) 

3xs(of)  3x5(a£)  3xs(o£)  ^  3x5(g£)  3x5(g£)  3x5(g£) 


8X^03)  3X2(o0)  3X3(0o)  j3x4(o0)  3xs(0o)  3xg(gg) 
3x6(of)  3x6(gf)  3x6(gf)  1 3x6(af)  3x6(af)  3x6(gf) 


BX^Oq)  3X2(o0)  3X3(03)  |  3x4(0q)  3x5(o0)  3x6(o0) 
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or 


SA  (oq) 

«x  (o0) J 


B-16x(of) 


These  equations  are  equivalent  to  Eq(Z-33)  and  EqQZ-34;  m  chapter  11. 

Numerical  Application  Of  The  Optimization  Algorithms 

The  application  of  the  numerical  optimization  techniques  of  Chap¬ 
ter  II  to  the  problem  derived  in  the  previous  section  is  shown  in  flow 
chart  form  in  Fig.  5,  on  the  following  page.  A  discussion  of  the  sol¬ 
ution  flow  is  presented,  with  the  considerations  given  for  the  choice 
of  each  algorithm. 

General  Discussion  Of  Algorithms.  To  start  the  solution,  the 
gradient  technique  is  employed.  The  determining  factor  for  this 
choice  is  that  since  the  gradient  is  not  defined  to  be  zero,  the  larg¬ 
est  convergence  envelope  is  obtained.  A  set  of  controls  are  chosen  to 
initialize  the  gradient  end  the  algorithm  proceeds  until  the  criterion 
for  convergence  is  satisfied.  At  this  point,  the  state  and  adjoint 
values  at  the  boundaries  are  stored  for  use  as  Initial  estimates  in 
the  terminal  error  function  technique. 

The  terminal  error  function  technique  is  used  as  an  intermediate 
technique  due  to  the  complexity  of  the  problem.  The  preliminary 
Investigation  into  other  problems  where  both  first  and  second  varia¬ 
tion  techniques  were  used  brought  out  one  major  common  problem.  The 
convergence  of  the  gradient  had  to  be  accurate  to  almost  ten  slgnfi- 
cant  figures  before  convergence  could  be  obtained  in  the  second  varia¬ 
tion  technique.  The  complexity  of  the  problems  surveyed  appeared  to 
be  less  than  the  complexity  of  the  problem  under  discussion  here, 
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therefore,  the  apparent  solution  waj  to  develop  an  intermediate  tech¬ 
nique  which  has  all  the  characteristics  of  the  second  variation  but 
none  of  the  perturbation  equations.  These  contain  the  higher  order 
effects,  and  are  the  most  sensitive  to  an  incompletely  converged  con¬ 
trol.  The  method  of  converging  to  the  boundary  conditions  improves 
the  rate  of  convergence  from  a  distance,  but  is  sluggish  when  the  de¬ 
sired  conditions  are  nearly  attained.  This  property  complements  the 
second  variation  convergence  property  very  well. 

The  boundary  conditions  on  one  end  of  the  problem  are  employed  as 
an  initial  guess  for  the  terminal  error  function,  these  are  refined  in 
an  iterative  fashion  until  the  second  variation  technique  accepts 
them.  The  complementary  convergence  property  is  employed  if  the  con¬ 
vergence  in  the  second  variation  appears  sluggish.  An  arbitrary  indi¬ 
cator  of  this  is  the  weighting  value  with  the  Newton-Raphson  matrix. 

The  second  variation  procedure  is  used  to  refine  further  the 
boundary  conditions.  The  procedure  will,  in  general,  converge  in  fewer 
iterations  than  the  other  techniques  since  it  is  a  second  order  method. 
However,  linear  assumptions  are  made  and  the  initial  estimates  to 
start  the  procedure  must  be  accurate  enough  to  prevent  divergence. 

When  the  method  converges,  the  convergence  is  generally  to  within 
eight  significant  digits  in  less  than  thirty  iterations. 

Some  numerical  considerations  unique  to  each  of  the  techniques 
are  discussed  in  the  following  sections. 

First  Variation  Numerical  Considerations.  During  the  thesis, 
there  occurred  several  areas  where  the  computational  solution  was  in¬ 
hibited  or  assisted  by  the  numerical  method  used.  Here,  some  of  the 
major  numerical  considerations  shall  be  discussed. 
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In  the  gradient  algorithm  there  exists  a  problem  in  the  integra¬ 
tion.  In  a  typical  optimal  control  problem,  there  are  two  integrations 
performed,  one  forward  for  the  states,  and  one  backwards  ror  the 
ad  joints .  Hers,  there  are  four  integrafionn.  These  are  necessitated 
because  of  the  sensitivity  coefficients  under  consideration,  and  the 
nature  of  an  entry.  The  functional  relationship  is  a  cumulative  one, 
such  that  all  sets  of  equations  functionally  related  to  the  one  being 
integrated  must  be  integrated  simultaneously.  The  table  below  indic¬ 
ates  the  four  distinct  sets  of  equations,  the  order  in  which  they  are 
Integrated,  and  the  location  of  the  boundary  conditions. 


Table  III 

Known  Boundary  Conditions  For  Gradient  Algorithm 


Integration 

Order 

Equation 

Set 

Known 

Boundaries 

1 

State 

Initial 

2 

Sensitivity 

Coefficients 

Final 

3 

Sensitivity 
Coefficient  Adjoint 

Initial 

4 

State  Adjoints 

Final 

The  method  of  successive  integration  sweeps  is  extremely  time 
consuming,  therefore  to  cut  down  a  large  amount  of  this  time,  the 
alpha  search  In  the  gradient  uses  only  the  first  two  Integration 
sweeps.  From  these  sweeps  the  cost  can  be  fully  determined,  and  the 
alpha  search  can  proceed  in  minimum  time.  There  is  a  pitfall,  which 
occurred  often, enough  to  warrent  a  check  in  the  computation  scheme. 
The  adjoint  equations  are  never  consulted  during  the  alpha  search, 
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and  sometimes  overflow  the  computer  due  to  the  large  changes  in  the 
cjntrol.  This  occurs  primarily  in  the  sensitivity  coefficient  adjoint 
equations,  which  are  third  order  equations,  when  cilia  uttuno,  the 
nearest  value  of  alpha,  for  which  a  gradient  could  be  calculated,  is 
used.  With  this  part  of  the  algorithm  streamlined,  a  great  amount  of 
computation  time  is  saved. 

In  the  terminal  error  function,  the  major  problem  area  is  in  the 
convergence.  This  problem  occurs  as  a  direct  res\ilt  of  the  functional 
relation  discussed  in  the  gradient  algorithm.  The  convergence  is  in 
general  slow  when  all  criterion  elements  are  balanced,  such  that  each 
boundary  condition  error  carries  a  similar  weight.  The  convergence  is 
much  improved  by  converging  on  each  element  In  turn,  saving  the 
adjoint  variables  until  last.  The  terminal  error  function  in  this 
problem  converges  very  quickly  to  the  state  boundaries,  at  the  expense 
of  the  adjoint  boundaries.  This  convergence  problem  will  be  further 
discussed  when  the  results  are  presented. 

Second  Variation  Numerical  Considerations.  In  using  the  second 
variation  method,  the  computer  time  required  becomes  of  primary  con¬ 
cern.  In  the  entry  problem  presented  in  this  report,  it  is  necessary 
to  solve  162  differential,  equations.  For  integration  of  a  trajectory 
with  100  Integration  intervals  using  a  fourth  order  Runge-Kutta  inte¬ 
gration  scheme,  approximately  one  and  one-half  minutes  of  computer 
execution  time  on  an  IBM  7040  is  required.  If  30  to  40  iterations 
with  the  Newton-Raphson  matrix  are  required,  a  total  execution  time  of 
about  45  to  60  minutes  results.  This  problem  Is  the  most  difficult  to 
overcome.  If  all  nine  position  sensitivity  coefficients  are  included 
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In  the  problem,  It  la  necessary  to  solve  24  state  and  adjoint  equations 
and  576  perterbation  equations  resulting  In  a  large  time  ir crease. 

Depending  on  the  nature  of  the  state  and  adjoint  equations,  the 
second  variation  procedure  may  converge  in  fewer  iterations  using  back¬ 
ward  integration  than  using  forward  integration.  However,  in  order  to 
reduce  the  terminal  velocity  by  several  thousand  meters  per  second  and 
the  terminal  flight  angle  by  siveral  degrees  while  holding  the  initial 
conditions  constant,  it  is  more  desirable  to  use  forward  integration. 
This  allows  the  Newton-Raphson  procedure  to  slowly  lower  the  terminal 
conditions. 


For  the  entry  problem  considered,  the  second  variation  procedure 
may  be  started  without  initial  estimates  from  either  the  gradient  or 
terminal  error  function  procedures.  This  requires  a  visual  examine- 
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er  signs  on  each  variable  estimated  to  cause  the  trajectory  to  proceed 


in  the  desired  direction.  The  magnitudes  can  be  determined  by  testing 


the  program  with  any  arbitrary  magnitudes  and  subsequently  reducing 


the  estimated  values  until  no  overflow  condition  exists  in  the  com¬ 


puter.  This  does  not  produce  an  acceptable  trajectory,  but  does  pro¬ 
duce  a  beginning  trajectory  to  allow  the  final  conditions  to  be  satis¬ 
fied  by  the  Newton-Raphson  procedure.  The  magnitudes  of  the  estimates 
for  one  example  case  are  presented  in  Appendix  B.  This  hueristic 
approach  is  not  mathematically  based  but  allows  maximum  use  of  engin¬ 
eering  judgement  to  start  the  procedure. 

The  Newton-Raphson  matrix  for  the  entry  problem  is  a  6x6  matrix 


and  must  be  inverted  with  double  presion  arithmetic  in  the  computer  to 
allow  convergence.  No  measures  of  ill-condition  for  this  matrix  have 
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been  computed,  but  the  elements  for  the  entry  problem  range  from 
10+l"  to  10-19.  It  is  possible  that  the  matrix  is  highly  ill-condit¬ 
ioned  because  of  this.  However,  on  some  of  the  trajectories  considered 
backward  integration  produced  convergence  to  within  10” ®  m  in  posi™ 
tion,  1O-0  mps  in  velocity,  and  1Q“5  deg  in  flight  angle. 
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V.  Trajectory  Solutions 

In  order  to  analyze  the  form  and  magnitude  of  the  sensitivity  co¬ 
efficients,  a  set  of  solutions  to  the  optimal  entry  problem  ia  prod¬ 
uced,  eacn  of  which  minimized  Lite  sensitivity  coefficients  for  a  given 
set  of  parameters  and  control  constraints.  In  this  chapter,  the  re¬ 
sults  of  these  solutions  are  presented  and  discussed. 

The  optimal  problem  is  cast  into  the  following  four  general  forms 

A.  Unconstrained  control  with  a  drag  control  bias 

B.  Unconstrained  control  without  a  drag  control  bias 

C.  Constrained  control  with  a  drag  control  bias 

D.  Constrained  control  without  a  drag  control  bias 
Within  each  of  these  forms,  the  boundary  conditions  are  varied  to  pro¬ 
duce  variations  in  the  lift  and  drag  controls  and  in  the  trajectories. 
The  bias  value  represents  the  results  of  defining  an  arbitrary  vehicle 
configuration.  The  concept  of  designing  the  controls  to  contain  all 
the  vehicle  parameters  is  described  in  Chapter  III.  The  bias  chosen 
is  the  mean  of  the  control  boundaries.  For  the  drag  control  element 
the  bias  value  is: 


( 


vi  ,  uDmax*’  UDmin 
%  L - o - 


0.00014688 


(5-1) 


The  lift  control  is  bounded  symmetrically  about  zero;  thus,  no 
bias  is  necessary. 

The  trajectory  solutions  are  presented  in  the  first  portion  of 
this  chapter,  and  are  discussed  in  the  latter  portion.  Seven  cases 
are  considered. 
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Unconstrained  Control  Solutions 

Four  trajectories  are  discussed  here,  each  with  a  different  con¬ 
trol  and  terminal  boundary  value  limitation.  The  first  two  cases  do 
not  include  the  arbitrary  configuration  bias  parameter.  The  second 
two  Include  a  defined  configuration  and  also  an  attempt  to  achieve 
nearly  horizontal  flight  with  a  terminal  velocity  of  about  Mach  5  at  a 
terminal  altitude  of  10  Km. 

Case  1.  Free  Terminal  Velocity  And  Flight  Path  Angle  And  A  Range 
Angle  Of  1.71  Degrees.  In  the  first  trajectory  considered,  the 
following  boundary  conditions  are  defined: 


Range  Angle 

Radius-  (r) 

Velocity  (v) 

Flight  Angle  (y) 

°0 

120  Km 

7420  mps 

120  degrees 

°f 

9  Km 

uncons trained 

unconstrained 

The  results  for  this  trajectory  are  shown  in  Fig.  6  through  Fig. 

9.  The  maximum  change  in  flight  angle  for  this  trajectory  is  0.13  deg. 
and  its  terminal  value  is  119.88  deg.  The  velocity  increases  monoton- 
ically  from  the  initial  value  to  a  terminal  value  of  7563  mps  (meters 
per  second).  The  reasons  for  the  nearly  constant  acceleration  are 
found  by  examining  the  values  of  the  lift  and  drag  controls.  As  the 
vehicle  enters  the  most  dense  portion  of  the  atmosphere,  the  controls 
begin  to  change  from  the  zero  level.  The  drag  control  becomes  more 
negative  so  as  to  cause  constant  acceleration  of  the  vehicle.  The 
trajectory  for  this  case  has  a  nearly  constant  flight  path  angle  indi¬ 
cating  that  the  velocity  sensitivity  coefficient,  [3r(<j£)/3v(o)  ] ,  is 
minimized  by  forcing  the  trajectory  to  be  as  near  a  straight  line  in 
space  as  possible.  For  a  fixed  Of,  the  control  for  this  type  of 
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trajectory  produces  smaller  sensitivity  coefficients  than  for  most  of 
trie  otner  cases.  However,  die  uuLui.i=  of  the  trajectory  limits  its 
practical  use  because  a  thrust  capability  is  required  to  produce  the 
negative  drag  profile. 

Case  2.  Specified  Terminal  States  And  A  Range  Angle  Of  1.71 
Degrees.  The  second  trajectory  considered  is  constrained  to  meet  ter¬ 
minal  boundary  conditions  on  all  three  states.  The  conditions  are 
as  follows : 


Range  Angle 

Radius  (r) 

Velocity  (v) 

Flight  Angle  (y) 

°0 

120  Km 

7420  mps 

120  degrees 

°f 

9  Km 

6500  mps 

119.88  degrees 

The  results  of  this  case  are  shown  In  Fig.  10  through  Fig.  13. 

Trtvffl  t’ru J  unt-ory  ra«ofn))1oo  t-Ho  ni*oy< a»jo  n  noo  < r»  f-Vi  r.r.jl® 

profiles  for  the  two  cases  are  similar.  The  effect  of  the  lower  vel¬ 
ocity  at  the  terminal  point  is  reflected  in  the  drag  control.  The 
drag  control  for  this  case  is  positive  throughout  the  trajectory.  . 
Above  the  altitude  of  12.3  Km,  the  control  is  below  the  range  of  most 
practical  non-thrusting  entry  body  controllers. 

In  this  trajectory  the  controls  do  not  deviate  from  zero  until 
the  dense  portion  of  the  atmosphere  is  entered,  at  approximately  35  Km, 
and  uearly  all  of  the  controlling  is  done  in  the  last  one-sixth  of  the 
trajectory.  This  is  reflected  in  the  velocity  profile,  which  remains 
nearly  constant  until  the  end  of  the  trajectory.  The  drag  control  is 
dominant  in  this  case  as  shown  by  the  lift  to  drag  ratio.  The  form 
and  magnitude  of  the  sensitivity  coefficients  is  similar  to  that  of 
Case  1. 
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The  maximum  deceleration  for  this  case  Is  138  G's  during  the 
final  portion  of  the  trajectory.  The  high  terminal  velocity  and  low 
drag  profile  In  this  case  limits  its  practical  use. 

Case  3.  Specified  Terminal  States  With  A  Control  Mean  And  A 
Range  Angle  Of  3.16  Degrees.  The  two  trajectories  discussed  previously 
are  limited  due  to  drag  control  values  below  the  boundaries  of  a  non¬ 
thrusting  controller.  To  eliminate  this  limitation,  two  additional 
elements  are  considered.  First  a  drag  control  bias  is  placed  in  the 
criterion  function  to  allow  for  drag  control  about  a  constant  value. 
Second i  the  terminal  conditions  of  the  problem  are  changed  such  that 
a  reasonably  flat  terminal  flight  angle  and  a  terminal  velocity  of  ap¬ 
proximately  Mach  5,  are  obtained. 

The  trajectory  for  this  case  has  the  following  boundary  conditions 


Range  Angle 

Radius  (r) 

Velocity  (v) 

Flight  Angle  (y) 

°0 

120  Km 

7  2 AO  mps 

108.0  degrees 

°f 

8.1  Km 

1650  mps 

98.9  degrees 

As  shown  in  Fig.  1A  through  Fig.  17,  the  controls  are  relatively  con¬ 
stant  until  the  vehicle  enters  the  most  dense  portion  of  the  atmos¬ 
phere  when  large  changes  begin  to  occur  in  the  velocity  and  y  profiles. 
The  controls  reach  maximum  values  at  the  maximum  dynamic  pressure 
point  which  is  at  10.5  Km.  The  lift  to  drag  ratio  is  also  a  maximum 
at  this  point.  The  changes  in  control  are  reflected  in  the  sensitiv¬ 
ity  coefficients,  especially  in  the  velocity  sensitivity  coefficient. 

A  secondary  objective  of  this  trajectory  is  to  reduce  the  amount 
of  deceleration  of  the  vehicle  by  obtaining  a  flatter  trajectory  and 
by  inducing  additional  drag  throughout  the  entry.  This  reduction  in 
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drag  deceleration  in  this  case  is  offset  by  the  lift  splice  shown  in 
the  graph  of  the  controls  resulting  in  a  high  total  deceleration  at 
the  end  of  the  trajectory.  Also,  during  the  highest  deceleration,  the 
drag  control  falls  below  acceptable  limits,  and  even  becomes  negative. 
Hie  negative  drag  requires  a  thrusting  vehicle,  and  thus  limits  its 
use.  Once  again  the  two  control  elements  do  not  change  significantly 
from  the  initial  values  until  the  last  portion  of  the  trajectory  in  the 
most  dense  part  of  the  atmosphere.  The  maximum  deceleration  is  101  G's 
which  occurs  at  approximately  the  point  of  maximum  dynamic  pressure. 

Case  4.  Specified  Terminal  States  With  A  Control  Mean  And  A 
Range  Angle  Of  3.32  Degrees.  In  the  previous  case,  the  sensitivity 
coefficients  appear  to  be  simply  a  function  of  the  range  angle.  To 
Investigate  this  relationship,  as  well  as  to  generate  a  flatter  tra¬ 
jectory,  an  additional  case  is  defined.  Hie  boundary  conditions  for 
this  case  are  defined  as  follows: 


Range  Angle 

Radius  (r) 

Velocity  (v) 

Flight  Angle  (y) 

120  Km 

6780  mps 

107.0  degrees 

°f 

8.1  Km 

1650  mps 

98.8  degrees 

The  results  of  this  case  are  shown  in  Fig.  18  through  Fig.  21.  The 
trajectory  obtained  is  similar  to  that  of  Case  3.  The  controls  for 
these  two  cases  differ  only  slightly  in  both  form  and  magnitude.  The 
sensitivity  of  radius  with  respect  to  velocity  for  this  case  rises 
sharply;  whereas,  the  sensitivity  of  the  radius  with  respect  to  flight 
angle  and  the  sensitivity  of  radius  with  respect  to  radius  appear  to 
continue  to  vary  linearly  with  range  angle  as  in  the  first  three 
cases. 
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For  the  secondary  objective  of  reducing  the  deceleration!  it  is 
shown  that  along  with  the  increase  in  one  sensitivity  coefficient  a 
lower  maximum  deceleration  is  attained.  The  maximum  deceleration  for 
this  trajectory  is  97.0  G's  occuring  again  at  approximately  the  point 
of  maximum  dynamic  pressure. 

Solutions  With  Lift  And  Drag  Control  Constraints 

All  of  the  four  previous  cases  satisfy  the  necessary  conditions  of 
the  optimal  entry  problem.  However,  none  satisfy  the  practical  re¬ 
quirement  that  the  controls  remain  within  physical  limits.  Even  in  the 
cases  with  a  drag  control  bias,  the  control  boundaries  were  violated. 

To  correct  for  this,  Valentine's  procedure  for  a  bounded  control 
solution  is  used  for  the  entry  problem  in  the  second  variation  opti¬ 
mization  method. 

Two  trajectories  are  determined  in  the  following  cases;  one  which 
employs  the  drag  control  bias,  and  one  which  does  not.  In  both  cases 
the  control  constraints  are  as  follows: 

2 .76*10"5  <1^5  2.20*10“3 

-1.6u_  5  u,  5  1.6u 

D  L  D 

The  lift  constraint  simply  limits  the  lift  to  drag  ratio. 

Case  5.  Specified  Terminal  States  And  A  Range  Angle  Of  2.0  Deg¬ 
rees.  This  trajectory  is  obtained  by  an  investigation  of  the  effects 
of  lower  starting  altitude,  velocity,  and  flight  angle  on  the  sensi¬ 
tivity  coefficients.  The  conditions  for  this  case  are  as  follows: 


Range  Angle 

Radius  (r) 

Velocity  (v) 

Flight  Angle  (y) 

°0 

115  Km 

6115  mps 

115.5  degrees 

0  Km 

5000  mps 

114.5  degrees 
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As  the  Fig.  22  through  Fig.  25  show,  this  trajectory  results  In 
large  changes  of  the  flight  angle  during  the  entry.  These  changes  are 
reflected  hv  control  curves,  which  show  that  nearly  all  the  maneu¬ 
vering  is  dona  by  the  lift  control ,  while  the  drag  control  remains  on 
the  lower  boundary.  The  lift  control  begins  to  significantly  deviate 
from  zero  at  about  the  midpoint  of  the  trajectory.  However,  this  does 
not  significantly  affect  the  velocity  and  flight  angle  until  the  last 
quarter  of  the  trajectory.  At  this  point,  lift  becomes  the  dominant 
control,  and  begins  to  directly  affect  the  flight  angle,  y.  Referring 
to  the  graph  of  the  lift  to  drag  ratio,  dominant  lift  control  does  not 
occur  until  the  lift  to  drag  ratio  approaches  unity.  The  lift  control 
reduces  the  flight  angle  to  well  below  the  terminal  value  and  the  '’rag 
control  increases  at  the  very  last  of  the  trajectory  to  satisfy  the 
terminal  velocity  condition  and  to  dive  to  the  terminal  flight  angle 
condition. 

The  drag  control  remains  on  the  lower  boundary  until  the  terminal 
phase  of  the  trajectory.  With  drag  remaining  at  a  minimum  for  nearly 
all  of  the  trajectory,  the  lift  control  maneuvers  the  flight  of  the 
vehicle . 

The  effect  of  the  control  constraint  is  evident  by  the  smooth 
changes  in  the  flight  angle.  The  lift  control  is  more  uniformly  dis¬ 
tributed  during  the  flight  than  in  the  previous  cases.  The  lift  and 
drag  controls  are  within  acceptable  practical  limits  throughout  the 
trajectory. 

This  trajectory  resembles  those  of  Case  3  and  Case  4  with  respect 
to  the  form  of  the  sensitivity  coefficients.  In  each  of  these,  the 
sensitivity  coefficients  appear  to  be  smooth  curves  until  changes  in 
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the  flight  angle  occur.  At  that  point  the  radial  and  velocity  sensi¬ 
tivity  coefficients  change  significantly.  This  effect  shall  be  dis- 

4*  kU*  pCttiCH  Z  f  thC  Cunptwlt 

This  trajectory  is  the  most  practical  yet  discussed.  However, 
due  to  the  high  terminal  velocity  and  the  high  deceleration  resulting 
from  the  high  drag  control  at  the  last  of  the  trajectory,  itB  use  is 
also  limited.  The  maximum  deceleration  is  83  G's  and  occurs  at  the 
terminal  point  of  the  trajectory. 

Case  6.  Specified  Terminal  States  With  A  Control  Mean  And  A 
Range  Angle  Of  2.06  Degrees.  This  case  contains  similar  boundary 
conditions  to  those  of  Cases  3  and  4,  but  with  a  smaller  total  range 
angle.  The  boundary  conditions  for  this  case  are  as  follows: 


Range  Angle  Radius  (r) I Velocity  (v)  Flight  Angle  (y) 


120  Km 
8.1  Km 


7667  mps 
1600  mps 


118.4  degrees 
97.0  degrees 


The  trajectory  for  this  case  is  very  similar  to  thoBe  found  in 
Cases  3  and  4.  The  results  are  illustrated  in  Fig.  26  through  Fig.  29. 
The  effect  of  the  control  constraint  is  evident  by  examining  the  por¬ 
tion  of  the  trajectory  during  which  the  flight  angle  changes.  The  use 
of  the  Valentine  procedure  keeps  the  controls  within  bounds;  however, 
this  alone  is  not  a  sufficient  practical  constraint.  The  rate  of 
change  of  the  lift  to  drag  ratio  caused  by  the  sudden  change  in  the 
magnitude  o<  the  drag  control  cannot  be  duplicated  in  any  practical 
vehicle. 

The  sensitivity  coefficients  generally  resemble  those  from  Cases 
3  and  4;  however,  there  are  some  notable  differences  between  these. 
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First,  due  to  changes  in  flight  angle  earlier  in  the  trajectory  than 
previously  noted,  the  radial  and  velocity  sensitivity  coefficients 
begin  tn  HptH ata  fro-  £  smooth  curve  curlier  is  the  trajectory  than  In 
Case  5.  Second,  they  appear  to  be  displaced  by  a  constant  value  from 
the  extrapolated  end  points  of  the  smooth  curves.  The  flight  angle 
sensitivity  coefficient  remains  a  nearly  linear  function  of  range 
angle.  There  is  also  in  each  Instance  where  a  flight  angle  change 
occurs  an  abrupt  change  In  slope  In  the  flight  angle  sensitivity  coe¬ 
fficient  curve.  In  Cases  3,  4,  and  6,  the  change  la  from  a  positive 
slope  to  a  leas  negative  slope,  while  in  Case  S  the  change  is  from  a 
positive  slope  to  a  more  negative  slope.  A  more  detailed  discussion  of 
these  tendencies  and  other  aspects  of  the  minimum  sensitivity  coeffic¬ 
ients  will  be  presented  in  the  last  section  of  this  chapter. 

This  trajectory  la  Impractical  due  to  rpto  of  oV>?r;*e  of  th® 
lift  to  drag  ratio.  The  maximum  deceleration  is  162  G's.  This  occurs 
at  the  point  of  maximum  control  change,  which  is  past  the  point  of 
maximum  dynamic  pressure.  It  is  evident  that  to  produce  a  practical 
trajectory,  it  Is  not  only  necessary  to  constrain  the  change  in  con¬ 
trol,  but  also  the  rate  of  change  of  control. 

Solution  With  Drag  Control  Constraints  Only 

It  has  been  previously  noted  that  when  drag  control  dominates  the 
entry  trajectory,  as  It  does  when  a  drag  control  bias  Is  used,  the 
changes  In  slope  of  the  sensitivity  coefficients  are  opposite  to  the 
changes  in  slope  when  a  bias  is  not  used.  These  changes  In  slope 
occur  only  when  the  flight  angle  changes  radically  during  the  trajec¬ 
tory.  Before  discussing  the  one  remaining  case,  It  should  be  noted 
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that  the  sensitivity  coefficient  equations  are  final  value  equations 
and  are  integrated  in  reverse  time.  For  this  reason  changes  in  the 
sensitivity  coefficients  are  more  easily  interpreted  by  considering 
them  beginning  at  the  final  range  angle  and  ending  with  the  initial 
range  angle.  Use  will  be  made  of  this  interpretation  method  in  the 
discussion  of  Case  7. 

Case  7.  Specified  Terminal  States.  Drag  Control  Constrained; 
Range  Angle  Of  1.86  Degrees.  The  results  of  the  trajectory  for  this 
case  are  illustrated  by  Fig.  30  through  Fig.  33.  The  boundary  condi¬ 
tions  are  as  follows: 


Range  Angle 

Radius  (r) 

Velocity  (v) 

Flight  Angle  (y) 

o0 

120  Km 

6700  mps 

120.3  degrees 

°f 

9  Km 

50C0  mps 

114.5  degrees 

The  trajectory  states  shown  in  the  figures  resemble  those  of  Case 
5,  with  two  significant  differences.  First,  the  amount  of  maximum 
change  in  the  flight  path  angle  is  twice  that  of  Case  5  due  to  the 
dominance  of  the  lift  control.  Second,  the  maximum  velocity  change  is 
increased  by  700  meter/second. 

The  lift  and  drag  profiles  of  the  two  cases  are  similar  if  magni¬ 
tudes  are  ignored.  In  this  case  the  lift  control  does  not  deviate 
significantly  zero  until  near  the  last  quarter  of  the  trajectory. 

The  drag  control  remains  on  the  lower  boundary  until  the  last  integra¬ 
tion  interval  of  the  trajectory.  Thus,  this  trajectory  consists 
predominately  of  lift  control. 

Hie  form  and  magnitude  of  the  radial  and  flight  angle  sensitivity 
coefficients  for  this  case  are  different  from  those  of  all  the  other 
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cases  considered.  The  radial  sensitivity  coefficient,  (3r(a^)/3r(a)], 
has  a  value  of  unity  at  the  terminal  point  and,  examining  the  coeffi¬ 
cient  backward  in  range  angle,  decreases  to  a  value  of  about  -0.05  and 
remains  at  this  value  until  the  initial  range  angle  of  zero  is  reached. 
The  flight  angle  coefficient,  [3r(o^)/3y(o) J,  starts  with  a  value  of 
zero  at  the  terminal  point,  decreases  to  about  -23000  and  then  increas¬ 
es  to  about  -9000  at  the  initial  point.  The  velocity  sensitivity  coe¬ 
fficient,  (3r(o£)/3v(o) ]  for  this  case  is  similar  in  form  to  those  of 
Case  5  .  The  coefficient  increases  from  zero  at  the  terminal  point  to 
a  value  of  about  +4.6  and  remains  at  this  value  until  the  initial  range 
angle  is  reached. 

This  trajectory  results  in  significantly  lower  range  and  flight 
angle  sensitivity  coefficients  than  for  any  of  the  other  cases.  The 
velocity  sensi fivi cy  coefficient  for  this  case  it  higher  then  fer  ell 
other  cases  except  Case  4. 

The  maximum  deceleration  for  this  trajectory  is  188  G's  at  the 
terminal  point  of  the  trajectory.  This  large  deceleration  is  primar¬ 
ily  due  to  the  change  in  control  during  the  dive  maneuver  during  the 
last  part  of  the  trajectory. 

Discussion  Of  Trajectory  Solutions 

In  the  previous  section,  seven  solutions  to  the  optimal  problem 
formulated  in  Chapter  IV  were  presented.  In  each  case,  different 
parameters  were  varied  to  allow  for  the  maximum  variation  among  the 
splutiona.  Although  the  solutions  varied,  there  were  some  general 
trends  which  were  noticeable  throughout.  In  this  section,  several 
observations  will  be  presented  and  substantiated.  The  observations 
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will  concern  the  curves  formed  by  the  Integration  of  the  sensitivity 
coefficients.  As  the  observations  shall  be  substantiated  by  examples 
from  the  cases,  a  summary  of  the  salient  points  from  each  case  arc 


Observation  1:  The  sensitivity  coefficient  solutions  may  be  sep¬ 
arated  into  two  distinct  portions  during  each  trajectory.  This  obser¬ 
vation  Is  substantiated  by  the  appearance  of  a  break  point  in  the 
three  sensitivity  coefficient  curves  in  all  cases  after  Case  2.  Ref¬ 
erring  to  Case  3,  there  is  a  smooth  curve  down  to  the  last  quarter  of 
the  trajectory,  then  the  curves  diverge  significantly.  The  first 
portion  of  the  trajectory  shall  be  defined  as  the  non-aerodynamic  por¬ 
tion.  The  second  part  shall  be  defined  as  the  aerodynamic  portion. 

Observation  2:  Hie  non-aerodynamic  portion  of  the  trajectory  la 
characterized  by  a  drag  domln?"'?*,  r*rd  small  changes  in  flight  angle. 
This  observation  is  substantiated  by  Cases  1  and  2  which  are  non-aero- 
dynamic  trajectories.  Both  of  these  exhibit  a  drag  dominance  in  the 
lift  to  drag  ratio,  and  have  a  small  change  in  flight  angle.  All 
other  cases  exhibit  this  same  tendency  throughout  the  non-  aerodynamic 
portion. 

Observation  3:  The  sensitivity  coefficient  curves,  in  the  non- 
aerodynamic  portion  are  either  functions  of  range  angle,;  os  ■constant. 
This  tendency  Is  common  throughout  the  cases.  This  observation  shall 
be  substantiated  for  each  sensitivity  coefficient.  The  ladial  sensi¬ 
tivity  coefficient  [3r(crf)/3r(a)  ] ,  in  each  case  is  a  near  constant  with 
respect  to  range  angle.  For  example,  in  Case  6,  the  value  Is  .7  from 
Oq  to  the  break  point.  The  velocity  sensitivity  coefficient, 
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[3r(o^)/3v(a) ] ,  In  all  cases,  excluding  Case  7,  is  observed  to  be  some 
exponential  function  of  range  angle.  In  Case  7,  the  value  appears  to 

•  .  i  f n _ ■ * J  i.t.  1>«  a 
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second  function  which  determines  the  magnitude  of  the  expontial  func¬ 
tion  of  range  angle.  This  function  remains  constant  throughout  the  ' 
non-aerodynamic  portion  of  the  trajectory.  The  flight  angle  sensi¬ 
tivity  coefficient,  [3r(crf)/3Y(o)  ] ,  is  a  linear  function  of  range 
angle.  In  all  cases,  except  Case  7,  the  slope  is  positive.  In  all 
cases,  however,  there  is  a  deflnate  linear  relation  to  range  angle. 

It  is  apparent,  from  the  cases  which  both  have  different  slope  change 
at  the  break  point,  that  there  is  some  function  which  sets  the  value 

of  — [3r(af)/3y(o) ]  equal  to  a  constant  at  the  break  point.  Since 
3o 

the  sensitivity  coefficients  are  integrated  in  reverse  time,  as  final 
value  equations,  the  conditions  for  tbs  ronstnnt  v»lui>  of 
(3r(a£>/3r(o) ] ,  the  constant  multiplier  in  3r/3v  and  the  first  deriva¬ 
tive  of  [3r(0f)/3y(o) ]  are  set  at  the  break  point  as  final  values  for 
the  reverse  time  integration.  These  conditions  shall  be  further  dis¬ 
cussed  after  the  observations  concerning  the  aerodynamic  portion  of 
the  trajectory  have  been  stated. 

Observation  4;  The  aerodynamic  portion  of  the  trajectory  is  char¬ 
acterized  by  changes  in  both  flight  angle  and  velocity;  and  dominance 
by  the  lift  control. 

This  observation  is  substantiated  by  the  fact  that  in  each  case 
which  has  two  distinct  portions  of  the  trajectory,  the  break  point  of 
all  of  the  sensitivity  coefficient  curves  occur  at  a  given  range  angle, 
at  that  range  angle,  in  the  profile  of  velocity  and  flight  angle  for 
the  same  case,  a  similar  break  point  can  be  found.  In  forward  time  at 
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the  break  point,  the  velocity  and  flight  angle  profiles  diverge  from 
their  previously  near  constant  state.  Also  at  the  range  angle,  for  the 
same  case,  the  lift  drag  ratio  passes  through  approximately  .:>.  This 
nhnwn  a  Hp.fMrvtf'fl  dnm/lnnnrfl  hv  lift  until  the  last  part  of  the  aero¬ 
dynamic  portion  of  the  trajectory.  * 

Observation  5;  In  the  aerodynamic  portion  of  the  trajectory,  the 
[3r(Of)/3r(o) ]  and  [3r(o£)/3v(o) ]  curves  are  a  function  of  the  rate  of 
change  of  velocity  and  Y«  The  (3r(aj)/3Y(o) ]  curve  is  a  function  of 
the  rate  of  change  of  lift  with  respect  to  o. 

The  first  part  of  this  observation  is  substantiated  by  correlating 
the  curves  for  [3r(0f)/3r(a) ]  and  [3r(a^)/3y(o) ]  with  the  velocity  and 
flight  angle  curves.  For  example,  in  Case  6,  when  [3y/3o]-0,  there  is 
also  a  zero  slope  achieved  by  the  [3r(0£)/3r(o)  ]  and  [3r(op/3v(o)  ] 
curves.  Also  in  this  same  case  there  is  s  pul  as  in  fh*»  velnMry 
at  a  range  angle  of  1.96,  and  a  corresponding  pulse  in  both  the  radial 
and  velocity  sensitivity  coefficient  curves.  Additional  examples  are 
found  by  correlating  the  curves  of  Cases  3  and  4. 

The  second  part  of  this  observation  is  substantiated  by  correlat¬ 
ing  the  lift  curve  slope  and  the  [3r(Oj)/3Y(o) ]  curve.  For  example, 
in  Case  7,  which  best  shows  the  correlation,  the  slope  of  the  flight 
angle  sensitivity  coefficient  curve  undergoes  a  change,  here  from  a 
positive  slope  to  a  negative  slope  at  the  same  range  angle  that  lift 
undergoes  a  slope  change.  In  other  cases  the  [3r(o£)/3Y(cr)  ]  slope 
change  is  not  as  dramatic  as  in  Case  7,  however  the  slope  change  can 
be  closely  correlated  in  a  similar  fashion  with  the  lift  slope. 

Observation  6:  There  exists  some  relation  between  the  aero¬ 
dynamic  portion  of  the  trajectories  and  the  conditions  during  the 
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non-aerodynamic  portion  lor  the  following: 

1.  The  constant  that  describes  the  radial  sensitivity  coef¬ 
ficient  i<srlOj;/ar(a)  J. 

2.  The  constant  which  multiplies  the  exponential  function 
that  describes  the  velocity  sensitivity  coefficient  [3r(o^)/Dv(o) ]. 

3.  The  slope  (sign  and  magnitude  of  the  first  derivative)  of 
the  linear  function  which  describes  the  flight  angle  sensitivity  coef¬ 
ficient  [3r(o£)/3Y(o) ]. 

Observation  7:  There  exists  some  relation  between  the  control 
boundaries  and  the  aerodynamic  portion  of  the  trajectory. 

Observation  6  and  7  are  existential  only  and  are  substantiated  by 
the  different  cases.  For  observation  6,  it  is  necessary  only  to  note 
the  variations  in  the  aerodynamic  portion  among  Cases  3  through  7,  with 
respect  to  the  non -aerodynamic  poiLlou  of  the  Uajeutoiy .  Iu  uu'uhluul- 
late  observation  7,  a  comparison  of  the  sensitivity  coefficients  for 
Cases  3  and  4  with  those  for  Cases  5,  6,  and  7  shows  some  relation 
exists.  No  further  extrapolations  can  be  made  for  observations  6  and 
7  due  to  a  lack  of  data  for  correlation. 

These  observations  are  intended  as  a  basis  for  conclusions  to  be 
drawn  concerning  the  general  problem  of  minimizing  the  sensitivity  co¬ 
efficients.  These  shall  not  be  drawn  here,  but  shall  be  drawn  follow¬ 
ing  the  implicit  guidance  simulation  presented  in  the  succeeding  chapt¬ 
ers. 
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VI .  Optimal  Implicit  Guidance  Simulation  For  The  Entry  Problem 

In  the  guidance  of  astronautical  vehicles,  there  are  two  basic 
philosophies  which  can  La  uscu>  Cue  plillusophy  is  called  explicit 
guidance  and  the  other  implicit  guidance.  These  guidance  philosophies 
are  described  in  considerable  detail  in  reference  2.  This  chapter  will 
define  implicit  guidance  and  optimal  implicit  guidance  and  will  outline 
the  procedures  for  simulating  the  optimal  implicit  guidance  philosophy 
on  a  high  speed  digital  computer.  The  purpose  of  such  a  simulation  is 
to  investigate  the  validity  and  limitations  of  the  solution  to  the 
optimal  open-loop  and  closed-loop  control  problems  in  a  practical  situ¬ 
ation  with  arbitrary  initial  errors  in  the  trajectory  states  and  white 
Gaussian  distributed,  zero  mean  observation  noise  present. 

Definition  Of  Optimal  Implicit  Guidance 

In  order  to  define  implicit  guidance,  It  Is  desirable  to  also  de¬ 
fine  explicit  guidance.  In  the  philosophy  of  explicit  guidance,  is  is 
assumed  that  an  approximate  closed  form  solution  of  the  guidance  equa¬ 
tion  exists  which  will  explicitly  relate  the  current  control  vector  to 
the  terminal  boundary  conditions.  Ihis  involves,  for  the  entry  prob¬ 
lem,  the  solution  of  a  set  of  transcendental  equations  at  each  control 
point.  A  Nowton-Raphson  algorithm  Is  generally  used  to  solve  these 
transcendental  equations.  However,  since  this  must  be  done  at  each  con¬ 
trol  point,  a  significant  amount  of  computer  time  is  necessary.  For 
this  reason,  another  scheme,  which  determines  the  control  vector  as  an 
implicit  fjnetion  of  the  terminal  boundary  conditions  is  generally 
used. 

In  implicit  guidance,  it  is  assumed  that  a  nominal  trajectory  and 
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control  are  computed  before  the  beginning  of  the  mission  (or  simulation 
in  this  case).  The  procedure  then  consists  of  applying  a  linear  theory 
to  errors  or  deviations  from  the  nominal  trajectory  in  order  to  satisfy 
the  terminal  boundary  conditions.  In  optimal  implicit  guidance  the 
nominal  trajectory  and  control  are  the  precomputed  open-loop  optimal 
trajectory  and  control.  The  linear  theory  in  optimal  implicit  guidance 
consists  of  applying  a  set  of  precomputed  linear  feedback  gains  to  the 
deviations  from  the  nominal  trajectory  so  that  a  neighboring  optimal 
path  is  followed  to  the  specified  terminal  boundary  conditions.  When 
the  actual  values  of  the  states  are  observed  with  observation  noise 


present,  a  Kalman  filter  (described  in  Chapter  II)  may  be  used  to 
obtain  optimal  estimates  of  the  actual  values  of  the  states  at  each 
control  point.  This  philosophy  is  illustrated  in  the  block  diagram  in 


Figure  34,  below.  In  the  diagram,  w  is  the  observation  noise  vector. 


Figure  34. 


Assumptions 


Optimal  Implicit  Guidance  Block  Diagram 


The  optimal  implicit  guidance  philosophy  is  based  on  the  following 
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assumptions: 

A,  The  deviations  of  the  actual  states  from  the  nominal 
states  are  within  the  linear  range  of  the  feedback  gains. 

R .  Hlfl  nnmlnnl  nrtrt  Cr/>1  nnit  fra  1  o/’ffl  tn  a  rrt  r/anffr — 

. .  wr  v-..-- .v-  — w  -  “J  *■/  N.w..w*« 

uous  and  not  adjacent  to  any  boundaries.  If  this  assumption  is  not 
made,  the  required  change  in  control  computed  from  the  feedback  gains 
could  cause  either  the  control  or  the  states  to  exceed  boundary 
limitations. 

C.  The  observation  noise  vector  consists  of  white  Gaussian 
noise  with  a  known  mean  and  covariance  matrix. 

D.  The  vehicle  dynamics  can  be  described  by  a  set  of  non¬ 
linear,  coupled,  first-order  differential  equations  with  deterministic 
coefficients. 


Entry  Feedback  Equations 

The  linear  feedback  scheme  uses  the  following  equation: 

«u(o)l 


6u(o)  ■ 


W(c) 


6x(o) 


(6-1) 


where  the  matrix  [  3u(cr)/  3x(o)  ]  is  derived  such  that  the  trajectory 
follows  a  neighboring  path  to  the  specified  terminal  conditions.  If 
all  six  states  of  the  entry  problem  are  considered,  then  the  matrix  is 
of  dimension  2x6  as  there  are  two  control  components.  However,  onboard 
a  practical  vehicle,  the  errors  in  the  sensitivity  coefficients  cannot 
be  computed.  For  this  reason,  the  feedback  scheme  uses  errors  in  the 
states  only.  It  is  assumed,  and  shall  be  demonstrated  that  with  no 
measurement  error,  the  feedback  scheme  will  satisfy  the  position, 
velocity,  and  flight  angle  terminal- conditions.  If  there  are  errors 
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in  measurement,  the  use  of  a  minimum  sensitivity  coefficient  nominal 
trajectory  will  aid  in  minimizing  these  errors.  This,  however,  shall 
not  be  demonstrated. 

The  equations  for  the  entry  problem  are: 


4un(o) 

3ud(o)  3ud(o)  3ud(o) 

8x^0)  8x2(0)  8x3(o) 

4uL(o) 

8ul(o)  8ul(o)  3ul(o) 

.8 Xj  (0 )  8x2(0)  8x3(0). 

4*1  (o) 

4xj  (o) 
4x3(o) 


(6-2) 


where  the  2x3  matrix  above  is  found  from: 


8u(o) 

M  — 

8J»(o) 

-1 

8G(a) 

+  3G(o) 

8^(o) 

3x(a) 

1 

-1 

3x(o) 

3_u(o) 

y 

8x(o) 

8X,(o) 

31  (of) 

3X(of) 

4 

and  where: 


3G(o) 

8u(o) 

1 

-1 

m 

l 

^UD 

n 

i 

0 

1 

Sr 

3G,(o) 

3G,  (0) 

3G,  (0) 

3G(o) 

am 

3xt(o  ) 

3x2(o) 

3x3(o) 

8x(0) 

3Gj  (0) 

3G2  (0) 

8Gj;  (0) 

3xj  (0 ) 

8x2(o) 

8X3(0) 

’  865 (0) 

3Gl(o) 

3Gt(0) 

8G(a) 

ma 

8Xi(o) 

8X2(o) 

3X3(0) 

31.(0) 

8G2(o) 

8G2(o) 

8G2(o) 

8Xl(o) 

3X2(o) 

8X3(o) 

(6-4) 


(6-5) 


(6-6) 
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‘  3 Mo) 

3*,(o) 

3*,(o) 

3Ho) 

3*1  (op 
3*2 (°) 

3*2  <°f) 

3*2 (0) 

3*3  (°p 
3*2(0) 

n 

3*iC°p 

3>3(c) 

"*2  <°f) 
3*3 (o) 

3*31° p 
3*3(0) 

3*1 (of) 

3*2  (op 

3*3  (op 

Sx^o) 

3xj(o) 

3xi (0) 

3x(o) 

3*i(of) 

3x2(o) 

3*2(°f) 

3x2 (o) 

3*  3  (op 

3x2 (0) 

3*  (op 

3*1 (of) 
3x3 (a) 

3*2(0  p 

3x3(o) 

3*  3  (°f) 
3x3(o) 

3*i(of) 

3*2  (Of) 

3*3  (Of) 

and  KUd  arr  the  weighting  factors  on  control  in  the  criterion  func¬ 
tion  and  the  elements  of  the  matrices  are  found  by  choosing  the  proper 
elements  from  the  s«1uHn«<»  to  pArfi.t+i  «**•■!««  <n  fhanrer 

IV  and  Appendix  A.  The  matrices  above  are  written  out  to  show  that 
they  are  not  the  6x6  matrices  from  the  perturbation  equations  but  only 
3*3  partitions  of  those  matrices. 


Digital  Computer  Simulation 

For  the  entry  problem,  two  simplifying  assumptions  are  made  in 
addition  to  those  made  in  the  previous  section.  These  are  as  follows: 

A.  No  system  (i.e. ,  plant  or  message)  noise  is  present. 
System  noise  could  be  considered  by  including  several  additional  terms 
in  the  filter  equations. 

B.  The  observation  vector  is  a  linear  function  of  the 
state  vector  and  contains  linear-additive  zero  mean  noise  with  a  con¬ 
stant  covariance  matrix. 
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The  purpoee  of  this  section  Is  to  present  the  algorithm  used  to 
simulate  the  optimal  Implicit  guidance  schema  on  a  high  spend  digital 
computer.  Fig.  35  shows  the  basic  indexing  method  used  for  simulating 
the  true  (i.e.,  nominal)  system  and  for  applying  the  Kalman  flltar. 


Figure  35.  Simulation  Indexing  Method 


In  order  to  simulate  the  guidance  scheme,  the  problem  must  first 
be  put  into  discrete  form  for  the  digital  computer.  In  the  algorithm, 
presented  below,  it  is  assumed  that  the  gain  matrices,  control  values 
and  nominal  states  for  each  point  of  the  trajectory  are  put  out  on 
binary  tape  by  the  program  which  generates  them. 

Step  1.  Based  on  a  priori  knowledge,  estimate  the  constant 
observation  noise  covariance  matrix,  W. 

Step  2.  Based  on  the  first  observation,  2^(i  *  0) ,  and 
some  arbitrary  scheme,  estimate  the  initial  optimal  estimate  ftg.  For 
linear  observations,  jeo  is  estimated  from 

Step  3.  Compute  the  initial  error  covariance  matrix,  To* 
For  linear  observations,  this  is 

r0  -  w 
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This  Is  computed  rather  than  estimated  since  the  equation  which  defines 
specifies  the  Initial  Kalman  gain  matrix,  KQ  (see  Chapter  II  on 
Initial  error  covariance). 

Step  4.  Using  the  precomputed  values  of  control  and  the 
known  initial  nominal  trajectory  conditio  is.  Integrate  the  state  equa¬ 
tions  forward  in  time  to  the  terminal  range  angle  to  define  a  new 
nominal  trajectory.  This  is  done  because  when  the  discrete  control 
values  are  used  as  forcing  functions  in  the  state  equations.  A  linear 
interpolation  must  be  used  in  the  Runge-Kutta  Integration  routine, 
resulting  in  a  slightly  different  nominal  trajectory  from  the  precom¬ 
puted  one. 

Step  5.  Simulate  the  initial  observation  with 

Xo  "  *orU*  +  Ho  (6-11) 

Where  Vq  is  Mn  **rn  in*m  ftatiaa-tan  nn'to*  wa^fni- , 

Step  6.  Integrate  tha  true  state  (from  xjrue)  to  the  (1+1) 
observation  point  using  a  step  else  h  by  solving: 


itrue  „ 

(6-12) 

where 

u  -  unom  +  6u 

(6-13) 

and 

6u  ■  [gain]  Sac 

(6-14) 

and 

Sx  -  xtrue  ~  xnom 

(6-15) 

It  Is  assumed  here  that  the  gain  at  any  point  will  be  associated  with 
the  complete  interval  immediately  following  that  point. 


Step  7.  Simulate  the  observation  noise  at  the  (1+1) 
observation  point  as: 

*i+i  -  *i+r + —i+i  <6-i6> 

Step  8.  Propogate  optimal  estimate  and  error  covariance 
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to  (1+1)  observation  point  to  form  the  best  estimates  3t(l+ljl)  and 
r(l+l  1)  given  everything  except  .  The  subscript  j  will  denote  the 
Intervals  between  1  and  (11*1).  A  step  size  of  H  is  used  for  these 
Suuiiitci. vals.  The  steps  in  this  procedure  are  as  followsi 

A.  .i  *  0 

B.  Define  x(j |j)  ■  and  r(j|j)  * 

C.  Solve  the  following  equations  to  obtain  the  beat 
estimates  at  the  (j+1)  subinterval: 


where 


A(j+l,j)  -  3lt~x'J(j"jt^  *  H  +  1 

(6-17) 

r(j+i|  j)  -  A(j+i.j)r<j|  j)a(j+iij)t 

(6-18) 

x,(J+l|j)-Itx(j+l|j),uj+1] 

(6-19) 

(6-70) 

u  j+1  »  4%  +  «JM+1 

(6-21) 

-  X(J+1|J)  -  4°n 

(6-22) 

D.  Go  back  to  C  with  J" J+1  noting  that  r(j| J)*T(J | J-l). 
Continue  the  solution  until  the  (11*1)  observation  point  is  reached. 

At  this  point: 

x(i+l|i)  -  x(j+l|j)  (6-23) 

and 

r(i+l|i)  -  r(J+l I J)  (6-24) 

Step  9.  Compute  the  following  quantities  at  the  (1+1) 
observation  point: 

A.  Kalman  Gain  Matrix: 


Ki+1  "  r(i+l|i)[r(i+l|i)  +  wf 


(6-25) 
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B.  Error  Covariance  Matrix: 

ri+l  ■  U-K1+1]r(l+l|l)[I-K1+1Jr  +  Ki+1r(l+l|i)Ki+1  <6- 

C.  Optimal  Estimate: 

*i+1  -  x(  1+1 1 1+1)  -  x(i+l|U  +  K1+1[xi+i  "  3c<±+ll±)  ]  (6- 

Step  10.  Compute  the  observation  error  and  estimation 
errors  with  respect  to  the  nominal  by: 

A*i+1  A  Zl+l  "  4+f 
A*i+1  A  ii+1  -  * j+1 

Step  11.  i  -  i+1 

Step  12.  Go  back  to  Step  6  and  continue  until  the  final 


range  angle  is  reached 
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VII .  Simulation  Results 

Simulation  results  for  two  of  the  open-loop  trajectories  consid¬ 
ered  in  Chapter  V  are.  selected  for  presentation  here.  These  are  for 
Case  1  and  Case  7.  Since  no  graphical  analysis  can  be  performed 
directly  on  the  estimated  states  because  of  the  magnitudes  of  the 
numbers  involved,  only  deviations  from  the  nominal  optimal  trajectory 
are  presented.  The  observation  and  estimation  error  quantities  are 
defined  sb  follows: 

4j £<o)  4y.(o)  -  x*(a)  -  [xtru4(o)  -  x*(a)]  +  w(a)  (7-1) 

and 

Ax(o)  4  x(®)  "  £*(<*)  (7-2) 

These  are  the  errors  in  the  observation  vector  and  in  the  optimal  esti¬ 
mate  with  respect  to  the  nominal  optimal  trajectory. 

For  each  case,  two  basic  sets  of  results  are  presented.  These 
are  as  follows: 

A.  Initial  trajectory  errors  of  +1000  m,  +50  mps,  and 
+0.2  deg.  (+0.4  deg.  for  Case  7)  in  r,  v,  and  y  respectively  with 
observation  noise  present. 

B.  Initial  trajectory  errors  of  -1000  m,  -50  mps,  and 
-0.2  deg.  (-0.4  deg.  for  Case  7)  in  r,  v,  and  y  respectively  with 
observation  noise  present. 

C.  Initial  trajectory  errors  of  +5000  m,  +100  mps,  and 
+1.0  deg.  in  r,  v,  and  y  respectively  both  with  and  without  observa¬ 
tion  noise  present. 

D.  Initial  trajectory  errors  of  -5000  m,  -100  mps,  and 
-1.0  deg.  in  r,  v,  and  y  respectively  with  observation  noise  present. 
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The  gain  matrix  elements  for  each  trajectory  are  presented  In 
Figures  36  and  37  and  lr.  Figures  63  and  64. 

The  unconstrained  control  results  for  the  cases  with  terminal 
velocities  of  1630  mps  and  1600  mps  are  not  considered  here  because  of 
large  amplitude,  high  frequency  variations  in  the  gain  matrix  elements 
during  the  final  portion  of  the  trajectory.  Due  to  the  discrete  na¬ 
ture  of  the  optimal  implicit  guidance  scheme,  a  smaller  integration 
step  size  would  be  necessary  to  produce  usable  gains  for  these  cases. 

It  should  be  noted  that  the  trajectory  of  Case  7  is  generated  with 
Valentine's  procedure  in  the  second  variation  program  for  a  drag  con¬ 
trol  constraint.  This  fact  violates  the  assumption  made  in  the  pre¬ 
vious  chapter  that  the  control  quantities  could  not  be  adjacent  to  any 
boundaries.  However,  since  this  trajectory  produces  the  lowest  sensi- 
Hv/ity  «.?  wrtht’ilo  to  examine  the  prapartica  in 

a  feedback  loop. 

Fig.  38  through  Fig.  62  show  the  errors  with  respect  to  the 
nominal  trajectory  and  the  error  covariance  matrix  (T)  elements  for 
Case  1.  Fig.  65  through  Fig.  89  show  the  corresponding  quantities  for 
Case  7.  In  the  graphs,  sigma  is  the  same  as  the  range  angle  and  has 
no  relation  to  the  standard  deviation  of  the  noise. 

In  all  examples,  the  graphs  show  that  the  feedback  gains  have  no 
effect  on  velocity  and  flight  angle  until  the  last  few  points  in  the 
trajectory.  In  Case  1,  100  gain  values  are  used  and  in  Case  7,  150 
gain  values  are  used.  Virtually  all  the  effect  of  the  gain  values  oc¬ 
curs  during  the  last  two-tenths  of  a  degree  of  range  angle.  If  the 
last  gain  value  is  not  used,  observations  have  shown  that  the  terminal 
values  are  in  error  by  as  much  or  more  than  the  initial  error  values. 
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Figure  38.  Case  1.1  -  Position  Error  of  1000  Meters 
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Figure  42.  Case  1.1  -  Off-Diagonal  Error  Covariance  Terms 
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Velocity  Error  of  -50  Meters/Second 
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Figure  49.  Case  1.3  —  Velocity  Error  of  100  Meters /Second  without  Boise 


Figure  50.  Case  1.3  -  Flight  Ang.e  Error  of  1.0  Degrees  without  Noise 
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Figure  51 .  Cose  1.3  -  Diagonal  Error  Covariance  Terms 
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Position  Error  of  5000  Metsrs 


Figure  .56.  Cane  1.4  -  Diagonal  Error  Covariance  Terms 


Figure  58.  Case  1.5  -  Position  Error  of  -5000  Meters 


Figure  59.  Case  1.5  -  Velocity  Error  of  -100  Meters /Second 


Figure  62.  Case  1.5  -  Off-Diagonal  Error  Covariance  Terms 


Figure  63.  Case  7  -  Drag  Control  Gains 


Flight  Angle  Error  of  -.4  Degrees 
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Figure  73.  Case  7.2  -  Diagonal  Error  Covariance  Terms 


Figure  74.  Case  7.3  -  Off-Diagonal  Error.  Covariance  Terms 


Figure  76.  Case  7.3  -  Velocity  Error  of  100  Meters/Second  without  Noise 


Figure  80.  Case  7.4  -  Position  Error  of  5000  Meters 
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Figure  89.  Case  7.5  -  Off-Diagonal  Error  Covariance  Terms 
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The  flight  angle  is  the  noat  sensitive  to  Initial  condition  errors 
from  the  nominal.  This  state  variable  is  brought  to  near  the  terminal 

ul^k  kti*  liaf  tun  o>jHn  vllilM.  Un  to  thil  n flint* .  the  flight 

angle  error  is  allowed  to  deviate  by  large  amounts  from  the  nominal. 

The  flight  angle  error  also  appears  to  be  the  most  sensitive  to 
the  addition  of  observation  noise.  This  is  shown  by  a  comparison  of 
the  errors  with  and  without  noise  present.  Overall,  the  flight  angle 
Is  most  sensitive  to  any  errors  in  state. 

The  position  and  velocity  errors  at  the  terminal  point  are  rela¬ 
tively  insensitive  to  the  initial. condition  errors  and  to  noise.  In 
all  the  examples  presented!  the  largest  error  in  terminal  position  is 
about  four  meters  and  the  largest  error  in  terminal  velocity  is  about 
six  meters  per  second.  Terminal  position  error  appears  to  be  almost 
completely  insensitive  to  the  magnitude  of  the  initial  errors. 

A  comparison  of  the  graphs  with  positive  and  negative  initial 
condition  errors  shows  that  the  linear  range  of  the  feedback  gains  is 
very  nearly  the  same  on  either  side  of  the  nominal  trajectory  with  only 
a  slightly  higher  peak  error  occurring  in  the  flight  angle  when  nega¬ 
tive  initial  condition  errors  are  present.  The  error  curves  with 
errors  on  both  sides  of  the  nominal  are  nearly  mirror  images. 

The  error  covariance  curves  for  the  examples  without  noise  present 
a.'e  shown  for  purposes  of  comparison.  The  same  initial  error  covar¬ 
iance  estimate  is  used  in  both  cases  along  with  the  same  noise  covar¬ 
iance  W.  When  no  noise  is  present,  the  off-diagonal  terms  of  the  error 
covariance  matrix  simply  reflect  the  coupling  between  the  state  equa¬ 
tions  and  the  validity  of  the  linearizing  assumptions.  A  comparison 
of  these  elements  with  and  without  noise  present  shows  that  the  off- 
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diagonal  terme  hava  vary  naarly  the  same  form.  Thia  indicate*  that 
the  off-diagonal  tarme  for  thia  aat  of  equations  has  very  little  rel- 

■  l<4nnaV4a  .bU  a  maJaa  M.  a..  —  _  —  — J.  ..  J..  e  .  i.«.  . 
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VIII.  Concilia  Iona 

The  basic  objectives  of  this  thesis,  to  investigate  a  eat  of 
optimal  solutions  and  to  produce  a  suboptlmal  feedback  scheme,  have 

haan  arroiiml  iahaM,  The  conclusion*  drawn  from  hVi*  specific  open—  end 

\ 

closed-loop  results  presentation,  as  well  as  conclusions  concerning  the 
general  problem  are  presented  in  this  chapter. 

Open-Loop  Problem 

The  following  conclusions  are  made  concerning  the  open-loop  re¬ 
sults  in  Chapter  V. 

A.  The  amount  of  lift  dominance  at  the  end  of  the  trajec¬ 
tory,  and  the  length  of  range  angle  to  go  when  the  lift  becomes  domin¬ 
ant  is  a  deciding  factor  In  minimizing  sensitivity  coefficients.  This 
concl  union  In  nnnnort-pH  hy  nhqprvaMnnn  5  and  fi  (Hnninnarl  nh  Hip  and 
of  Chapter  V. 

B.  The  appearance  of  minimum  control  for  the  first  portion 
of  the  trajectory,  stems  from  the  presence  of  control  in  the  criterion 
to  be  minimized. 

Closed-Loop  Problem 

The  following  conclusions  are  made  concerning  the  closed- loop 
results  in  Chapter  VIZ. 

C.  The  implicit  guidance  law  simulation  indicates  that 
control  for  the  trajectories  tested  can  be  maintained  even  with  large 
initial  condition  errors,  and  will  drive  the  vehicle  along  a  neigh¬ 
boring  trajectory  to  the  boundary  conditions  with  little  error.  In 
this  simulation,  the  flight  angle  seems  to  be  the  most  sensitive  to 
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initial  errors  in  tha  states.  For  cha  caaaa  In  which  tha  velocity 
dropa  to  laaa  than  1700  ops  at  tha  terminal  point,  tha  feedback  gaina 
contain  such  high  frequency  cun  cent  during  cne  iaat  of  tna  flight,  that 
the  sampling  rate  is  not  high  enough  to  allow  the  gaina  to  be  need: 

0.  Tho  position  error  is  least  susceptible  to  both  pertur¬ 
bation,  and  observation  noise.  This  is  attributable  to  the  fact  that 
the  nominal  trajectory  is  a  minimum  position  sensitivity  coefficient 
trajectory.  This  confirms  the  concept  of  minimizing  sensitivity  coef¬ 
ficients  to  achieve  minimum  terminal  position  error. 

General  Problem 

The  following  conclusions  are  made  concerning  the  overall  problem: 

E.  The  use  of  a  passive,  non-thrusting  control  in  a 
classic  optimal  control  problem  is  impractical  unless  it  can  be  removed 
from  the  criterion  function.  In  this  problem,  the  cost  of  control  was 
a  deciding  factor  for  the  optimality  of  the  trajectories. 

F.  The  use  of  lliu  and  drag  for  control  requires  a  dense 
atmosphere,  and  thus  a  flat  trajectory  in  the  dense  region,  to  achieve 
a  well  controlled  entry. 

G.  The  minimum  sensitivity  coefficient  trajectory  is  one 
in  which  the  vehicle  begins  a  nominal  entry  angle,  levels  out  to  nearly 
horizontal  flight,  and  then,  during  the  last  part  of  the  trajectory, 
dives  to  meet  the  terminal  angle  condition.  This  appears  to  confirm 
conclusion  F,  and  is  illustrated  by  Case  7. 

H.  For  each  of  the  trajectories  additional  constraints, 
other  than  those  contained  here  in,  should  be  employed  to  make  the 
trajectory  usable.  This  is  due,  in  part,  to  the  .type  of  control  used, 
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uhleh  works  only  during,  at  bast,  tha  last  quartar  of  tha  trajectory. 
For  soma  trajectories,  additional  elements  are  necessary  to  achieve  the 
magnitude  and  rate  of  change  of  magnitude  necessary  m  tne  optimal 
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IX.  Recommendations 

_ .  •  .  a  _  i lI -  — .  _  —  1  .««L1  AM  n  C  mt ina4  fllll  fV 
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coefficient  minimization  for  an  entry  trajectory  and  suggestions  for 
areas  of  further  study  are  presented  in  this  chapter. 

Sensitivity  Coefficient  Minimization 

The  general  problem  of  sensitivity  coefficient  minimization  for  an 
entry  trajectory  is  highly  dependent  on  the  lift-to-drag  ratio  and  the 
length  of  the  trajectory  during  which  control  is  exerted.  The  follow¬ 
ing  recommendations  are  made  in  this  area: 

A.  This  problem  should  be  eliminated  aB  a  problem  in 
classical  optimal  control  theory  and  cast  as  a  parameter  sweep  problem. 
The  purpose  of  this  is  to  sweep  only  within  a  range  of  practical  entry 
controllers . 

B.  The  parameters  to  be  considered  in  the  sweep  Include: 

1.  Lift-to-drag  ratio 

2.  Rate  of  change  of  lift-to-drag  ratio 

3.  The  range  angle  arc  for  which  the  lift  control  is 
allowed  to  dominate  the  trajectory. 

Suggestions  For  Areas  Of  Further  Study 

In  light  of  the  recommendations  given  above,  the  following  areas 
are  suggested  for  further  study: 

A.  For  the  parameter  sweep  problem: 

1.  Use  numerical  analysis  to  describe,  in  a  least 


square  sense,  smooth  functions  identified  during 
the  non-aerodynamic  portion  of  the  trajectory. 
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2.  Also  use  numerical  techniques  to  attempt  to  iden¬ 
tify  functions  which  may  be  used  to  describe  the 
changes  in  the  sensitivity  coefficients  during 
the  aerodynamic  portion  of  the  trajectory. 

3.  Identify,  using  the  relationships  determined  in  1 
and  2  above,  a  relationship  which  describes  the 
sensitivity  coefficients  throughout  the  trajec¬ 
tory. 

4.  Investigate  the  use  of  a  hybrid  computer  for  the 

parameter  sweep  problem  and  determine  its  value  in 
solving  an  optimal  control  problem  of  this  type 
with  the  primary  objective  of  obtaining  a  prac¬ 
tical  controller.  NOT  REPRODUCIBLE 

•  i 

Vrnnnfl  f*Vi a  nraa art f*  nmUImm  Tnnl  •'f 

•  a  *  .......  —  —  —  —  -  — —  **"*•"’  *-*- **  ** 

freedom  to  determine  the  effects  of  cross-range  motion  on  the  sensi¬ 
tivity  coefficients. 

C.  Investigate  the  present  problem  with  constraints  incl¬ 
uded  on  the  following  variables  to  achieve  a  practical  controller: 

1.  Drag  control  for  deceleration  constraints 

2.  Lift-to-drag  ratio  for  configuration  constraints 

3.  Rate  of  change  of  drag  control  for  maneuvering 
capability  constraints. 

4.  Rate  of  change  of  lift-to-drag  ratio  for  maneu¬ 
vering  capability  constraints. 

5.  Velocity  for  heating  constraints 

6.  Rate  of  change  of  velocity  for  deceleration 
constraints. 
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7.  Rate  of  change  of  flight  angle  for  deceleration 
constraints. 

uuti.  uic  uptiouil  imjfilitii,  {juiuunce  acheme  oc  iaveuti- 

gated  to  determine  the  most  practical  method  of  Implementing  it  and  the 
simplest  suboptimal  estimation  technique  which  can  be  used  in  an  actual 
entry  trajectory. 
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Appendix  A 

Open-Loop  Optimal  Control  Equations  Tor  The  Entry  Problem 
The  purpose  of  this  appendix  is  to  present  the  detailed  equation* 
fur  both  first  and  second  variations  of  the  optimal  entry  problem. 

They  are  presented  in  order  of  development. 
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Appendix  B 


Parametera  for  a  Sample  Trajectory 


The  purpose  of  this  appendix  is  to  present  some  sample  parametera 
to  be  used  In  starting  a  trajectory  solution.  The  'general  criterion 
function  for  the  entry  problem  is  (no  control  biac  included): 


For  one  particular  case,  that  of  Case  7  presented  in  Chapter  V,  the 
parametera  after  convergence  are  as  follows: 

K4  -  1.0  x  io"iO 

K5  -  1.0  x  10"10 

Kg  -  1.0  *  10"*° 

KUJ)  -  1.0  x  103 

K  -  1.0  *  103 

UL 

Mo.  Intervals  -  150 

^(oq)  -  -0.21064992  xio"11  M°f>  “  -0.53273489  x  10““' 

X2(o0)  -  -0.21462654  xio"10  Xz(of)  -  40.26407696  x  l0_n 

X3(o0)  -  40.59012141  xl0~°6  X3(af)  -  -0.56110516  x  10’07 

X4(of)  -  -0.59956528  x  lo'*0 

\5(of)  -  -0.54781552  x  io~“ 

\c(o e)  -  0.83000869  x  l(T15 
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